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Preface 


The  success  of  the  Space  Shuttle  Transportation  System 
will  have  a  great  deal  of  influence  on  this  nation's  policy 
and  attitude  toward  space.  Since  the  U.S.  Air  Force  has  been 
given  charge  of  military  operations  of  the  Space  Shuttle, 
opportunities  in  terms  of  operations  and  support  for  this 
versatile  system  are  increasing  in  the  military,  I  am 
personally  very  interested  in  being  a  part  of  this  new  mission 
of  the  Air  Force.  Because  of  my  great  personal  interest  smd 
desire  to  learn  more  of  the  Shuttle's  capabilities,  I  wel¬ 
comed  any  opportunity  to  involve  my  thesis  work  in  a  related 
subject,  with  an  adviser  well  informed  on  space  Shuttle 
operations,  I  was  extremely  fortunate  when  Capt  James  K, 

Hodge  was  assigned  to  the  faculty  here  at  AFIT.  With  him 
came  the  chance  to  learn  much  more  of  the  Space  Shuttle 
system  and  of  flight  into  and  from  the  space  environment. 

The  completion  of  this  project  would  not  have  been 
possible  without  the  guidance,  support,  and  tolerance  of  a 
number  of  people.  The  effort  I  put  forth  would  have  been 
futile  without  the  constant  guidance  smd  supervision  of  Capt 
Hodge,  who,  as  my  adviser,  steered  me  through  many  obstacles 
and  delays,  smd  led  me  to  the  results  found  here.  In  addi¬ 
tion,  Dr.  Wilbur  Hankey  was  very  helpful  in  directing  me  to 
resources  smd  in  raising  questions  that  could  only  be  posed 
through  experience,  I  wsmt  to  also  thank  Dr,  Jsimee  Hitchcock 
for  being  available  to  me  at  all  times  and  providing  addition- 
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t 
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al  methods  to  me*  Dr.  Joe  Shang  and  Mr,  George  Havener  also 
deserve  thanks  for  the  computer  program  used  in  this  investi¬ 
gation;  Lastly,  I  want  to  especially  thank  my  wife  Lisette 
for  patience  and  support  as  I  fretted  my  way  through  this 
project,  and  for  her  help  in  preparing  this  report, 

Peter  5.  Cappelano 
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Abstract 

V/iiid  tunnel  test  runs  v:ere  conducted  on  Space  shuttle 
orbiter  insulating  articles.  The  data  from  one  HRSI  run  was 
processed  by  a  heat  estimation  program,  revealing  a  large 
discrepancy  in  magnitude  of  the  heat  transfer  coefficient  over 
the  HRSI  tile  v;hen  compared  to  that  yielded  by  flat  plate 
theory  and  thin  skin  results.  An  initial  suggestion  for  the 
cause  of  this  discrepancy  was  the  nonisothermal  wall  effect. 
This  effect  was  the  subject  of  this  investigation.  The 
boundary  layer  heat  transfer  coefficient  v^as  evaluated  using 
a  numerical  boundary  layer  program,  supplemented  by  wall  tem¬ 
perature  superposition  analysis.  The  program  v:as  shovm  to 
give  results  consistent  with  Eckert  high  speed,  flat  plate 
theory  for  the  isothermal  wall.  A  comparison  of  Sutherland's 
Lav;  and  a  lov;  temperature  viscosity  calculation  showed  little 
difference  for  the  high  speed  flov;.  Theory  based  on  reference 
temperature  supported  this  result.  The  program  was  modified 
for  wall  temperature  step  over  the  insulating  article  to  model 
heat  transfer  coefficient  dependence  on  wall  temperature  vari¬ 
ation.  The  step  size  v;as  varied  by  doubling  and  quadrupling 
the  initial  step  size,  and  by  doubling  the  initial  wall  tem¬ 
perature.  The  magnitude  of  the  heat  transfer  coefficient  de¬ 
creased,  and  response  slowed,  with  increasing  step  size.  The 
step  location  v;as  also  varied.  Analysis  by  wall  temperature 
superposition  yielded  results  similar  to  program  results  with 
the  initial  wall  temperature  step.  The  nonisothermal  wall 
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effect  appears  to  account  for  most  of  the  discrepancy  in  the 
original  wind  tunnel  test  results. 
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HEAT  TnANS?':3^/B0UNDARY  LAYER 
:  INVESTIGATION  OP  HEATING  DISCREPANCIES 

■  IN  V/IND  TUIHTEL  TESTING  OP  ORBITER  INSULATING 

ARTICLES 

I.  Introduction 

Background 

The  successful  introduction  of  the  Space  Shuttle  system 
into  service  marks  the  beginning  of  a  new  era  in  the  utiliza¬ 
tion  of  the  space  environment,  A  manned  vehicle  capable  of 
entry  into  space  and  controlled,  gliding  reentry  with  the 
additional  quality  of  reusability  is  a  great  step  toward  in¬ 
creased  affordability  and  reliability  in  terms  of  space  trans¬ 
portation.  Repeated  use  translates  directly  into  reduced  cost 
to  travel  into  space.  Controlled,  gliding  reentry  means  greater 
versatility  in  selection  of  landing  sites  and  reentry  points. 
These  two  qualities  separate  the  space  Transportation  System 
(STS)  from  previous  modes  of  space  travel.  These  qualities 
also  place  greater  emphasis  on  the  capability  of  the  Qrbiter 
Thermal  Protection 'System  (TPS)  in  the  reentry  heating  en¬ 
vironment  (Ref  1  :2).  Reusability  requires  that  the  TPS  be 
able  to  sustain  reentry  heating  repeatedly.  Gliding  entry 
requires  that  the  TP3  be  able  to  sustain  reentry  heating 
through  a  variety  of  maneuvers  and  attitudes.  Consequently, 
the  capability  of  the  Thermal  Protection  System  must  be  thor¬ 
oughly  understood  through  testing  and  evaluation. 


The  Air  Force  Plight  Test  Center  (APFTC)  has  been  con¬ 
ducting  this  evaluation  through  aerothermodynamic  performance 
simulation,  and  collection  and  reduction  of  in-flight  thermo¬ 
couple  data  from  the  Reusable  Surface  insulation  (RSI)  on  the 
Orbiter.  The  data  reduction  was  accomplished  using  a  computer 
code  to  extract  aerodynamic  heat  rate  and  TPS  parameters  from 
transient  flight  test  maneuvers.  This  program  is  referred  to 
as  H3AT  ESTimation  or  HEATEST.  In  addition  to  data  from  ac¬ 
tual  Shuttle  flights,  thermocouple  measurements  from  Mach  14 
wind  tunnel  tests  v;ith  TPS  articles  were  obtained  and  re¬ 
duced  (Ref  1:2), 


Problem  Statement 

The  v.'ind  tunnel  testing  v/as  conducted  by  the  Air  Force 
Rright  Aeronautical  Laboratories  (APRAL)  in  a  circular  cross- 
section,  Mach  14  v;ind  tunnel.  Tests  vrere  made  using  three 
different  type  articles:  a  thin  skin  stainless  steel  plate, 
a  High-temperature  Reusable  Surface  insulation  tile  (HR3I1, 
and  a  Flexible  Reusable  Surface  Insulation  tile  (FRSIl.  The 
test  articles  v/ere  flush  mounted  on  a  test  plate  with  a  wedge- 
shaped  leading  edge.  'I'hermocouple  measurements  were  taken  at 
several  streamv;ise  locations  for  each  article,  as  indicated 
in  Figs  1,  2,  and  3*  A  total  of  59  test  runs  v;ere  catalogued. 
Runs  1-26  were  conducted  with  an  1R3I  tile.  During  these  runs, 
the  deflection  angle  of  the  test  plate  was  varied.  This  served 
tvio  purposes.  First,  the  angle  could  be  made  large  to  induce 
faster  heating  response  and  subsequently  lowered  in  order  to 
achieve  heating  equilibrium  on  the  test  article.  This  would 


reveal  a  limiting  temperature  for  given  conditions.  Second, 
variation  of  the  deflection  angle  could  simulate  flight 
maneuvers  and  give  some  indication  of  transient  heating  re¬ 
sponse,  The  actual  sequences  used  are  given  in  Table  I,  In 
addition.  Runs  27-45  were  conducted  with  an  PRSI  test  article 
and  Runs  46-59  were  conducted  vrith  a  thin  skin  article.  The 
sequences  for  these  are  given  in  Tables  TI  and  III,  respec¬ 
tively,  As  indicated  in  some  cases,  a  shock  generator  was 
used  to  simulate  shock  interaction.  Also,  in  the  HRSI  and 
FRSI  test  i^ns,  a  water  cooled  plate  parallel  to  the  test 
articles  was  used  to  maintain  a  somewhat  reliable  radiation 
sink  temperature  (Ref  1:8), 

’.7ind  tunnel  data  for  test  run  4,  with  an  KRSI  test  arti¬ 
cle,  was  input  to  the  analysis  program,  flEATSST,  to  extract 
heat  transfer  data.  The  sequence  of  variation  of  the  deflec¬ 
tion  angle  for  this  run  and  the  associated  temperature-time 
history  are  given  in  Pig  4,  The  plate  v;as  injected  at 
0L=  14  degrees  and  held  for  approximately  ten  seconds,  allow¬ 
ing  rapid  temperature  response.  The  deflection  angle  was  then 
lov:ered  to  oC=  5  degrees  and  the  temperature  appeared  to 
level  out,  indicating  that  thermal  equilibrium  was  achieved 
at  this  deflection  angle.  Finally,  the  plate  was  returned  to 
OC  =  14  degrees  and  the  temperature  climbed  for  the  remainder 
of  the  run.  The  test  indicated  that  equilibrium  was  achieved 
at  (3^=  3  degrees  . 

The  data  analysis  accomplished  using  H3ATRST  is  given  in 
Fig  5«  In  this  figure,  the  indicated  results  reduced  through 
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HSATTilST  differ  markedly  from  Eckert  flat  plate  theory  and 
thin  ^in  results.  At  the  deflection  angle  flC=  3  degrees  , 
theory  indicates  that  the  expected  value  of  the  heating  ratio 
should  be  =  “'•7  •  '*’he  KEATTDST  data  indicates  the 

value  to  be  h/h^^^  =0.95  •  Additionally,  an  equilibrium 
heat  rate  calculation  assuming  a  radiation  sink  temperature 
of  T  =  0  R  agrees  with  the  HEAT33T  value.  The  sane  cal¬ 
culation  assuming  a  sink  temperature  of  T  =  533  R  gives  a 
value  of  h/h^g^  lower  than  that  derived  from  HEATE3T  ^Ref  1:9). 

Prom  this  data,  the  reliability  of  the  wind  tunnel  tests 
is  questioned.  A  discrepancy  is  seen  between  HR3I  data  re¬ 
sults  and  Eckert  theory.  Thin  skin  test  data  seem  to  agree 
with  theory,  indicating  an  effect  only  in  the  IIRSI  case. 
Therefore,  the  primary  causes  indicated  are  those  differences 
betv;een  thin  skin  and  liRSI  materials.  The  ’JRST  material  has 
thermal  properties  quite  different  from  the  thin  skin  stain¬ 
less  steel.  Due  to  extremely  low  conductivity,  the  HT.SI  tile 
surface  temperature  should  increase  much  more  rapidly  than  the 
stainless  steel.  At  a  given  angle  and  time,  the  rasi  tile 
should  be  much  hotter  than  the  stainless  steel  surface  up¬ 
stream,  The  result  is  a  discontinuity  in  temperature  from 
the  steel  plate  to  the  1IHSI  tile.  This  is  the  nonisothermal 
v;all  problem.  Additionally,  the  HRSI  tile  achieves  equili¬ 
brium  because  it  has  hi_;h  emiasivity.  The  stainless  steel 
does  not,  A  valid  estimate  of  the  radiation  sink  temperature 
is  essential  to  the  accuracy  of  the  equilibrixim  heating  cal¬ 
culation.  Although  a  parallel  plate  of  known  temperature  is 
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used  as  a  sink,  the  tunnel  wall  temperature  is  unknown  and 
must  be  considered  as  a  sink  for  radiation,  based  on  view 
factor;  Additional  causes  of  error  that  must  be  considered 
include  inaccuracy  in  data  reduction,  inaccuracy  in  thermo¬ 
couple  data,  and  variation  in  the  boundary  layer  and  surface 
heating  ahead  of  the  test  articles  where  no  data  was  taken. 

Objectives  and  Aa-oroach 

In  order  to  make  the  v^ind  tunnel  test  data  valid  and  usa¬ 
ble,  it  must  be  corrected  to  account  for  the  possible  inac¬ 
curacies  cited.  The  greatest  effects  expected  are  those  due 
to  the  nonisothermal  v;all  problem,  and  those  caused  by  mis¬ 
understood  radiation  effects,  'fithout  knovrledge  of  the  tunnel 
wall  temperature,  the  radiation  cannot  be  accurately  estimated. 
This  precludes  pursuing  a  thorough  investigation  of  thermal 
radiation  and,  by  equilibrium  analysis,  derivation  of  heating 
data  on  the  test  article.  By  necessity,  the  focus  of  this 
investigation  v/ill  be  the  nonisothermal  wall  effect,  and  the 
resultant  change  in  heating  attributed  to  it.  Through  this 
investigation,  the  portion  of  the  discrepancy  due  to  the  non¬ 
isothermal  v/all  alone  may  be  determined.  Additionally,  at¬ 
tributing  the  discrepancy  cited  to  causes  unique  to  the  wind 
tunnel  should  render  greater  validity  to  theoretical  estima¬ 
tion  of  HRSI  heating.  Knowing  the  effect  of  the  nonisother¬ 
mal  wall  would  allov;  for  some  estimation  of  the  combined  ef¬ 
fects  of  the  other  sources  cited,  Finally,  the  effects  of 
the  nonisothermal  wall  condition  can  be  related  to  numerous 
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applications,  including  the  Space  Shuttle  ortiter.  several 
applications  v;ill  be  discussed. 

Analysis  of  the  nonisothermal  v^all  effect  will  require 
understanding  of  the  boundary  layer  development,  and  the  ef¬ 
fects  of  heating  on  boundary  layer  grov;th.  This  analysis 
will  best  be  accomplished  through  boundary  layer  solution 
techniques.  The  three  methods  of  solution  considered  here  are 
(1^  boundary  layer  similarity,  (2)  solution  of  the  boundary 
layer  equations  by  finite-difference  technique  in  a  boundary 
layer  program,  and  (3)  wall  temperature  superposition  analysis 
The  similarity  solution  of  the  boundary  layer  was  found  to 
be  invalid,  as  discussed  in  Appendix  A.  The  results  of  the 
boundary  layer  program  and  wall  tenpera.ture  superposition 
method  are  presented  subsequently.  The  boundary  layer  pro¬ 
gram  is  discussed  in  Appendix  B. 


II.  Theory 

Sckert-plat  Plate  Theory 

Heat  transfer  over  a  flat  plate  can  be  calculated  simply 
by  Eckert  theory  based  on  several  assumptions.  It  is  assumed 
that  the  flov^  medium  is  an  ideal,  single  component  gas,  with 
no  dissociation.  The  flow  over  the  plate  is  assumed  to  have 
locally  uniform  velocity  and  direction.  The  boundary  layer 
is  assumed  thin  enough  that  pressure  is  transmitted  without 
change  to  the  surface,  surface  temperature  is  assumed  con¬ 
stant,  Also,  the  fluid  is  assumed  to  have  constant  proper¬ 
ties,  independent  of  temperature  and  pressure.  The  following 
dimensionless  quantities  are  defined  as 

Nu^  =  h^  X  (1) 

St^  =  _!!x _  (2) 

^  P  c  u 
re  p  e 

where  h^is  the  heat  transfer  coefficient,  k  is  conductivity, 
p  is  density  of  the  fluid,  c^  is  the  specific  heat  at  con¬ 
stant  pressure,  and  u  is  the  velocity  of  the  fluid  at  the 
edge  of  the  boundary  layer,  ITu^  is  the  local  NUsselt  number 
and  St^  is  the  local  Stanton  Number.  For  laminar  flow  over 
the  plate,  the  following  relations  define  the  temperature 
recovery  factor  r,  the  Nusselt  Number,  and  the  Stanton  Num¬ 
ber: 

r  =/pr  (5) 
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(4) 


NU, 


= 


Ku, 


Pr  Re, 


=  0.332  Pr’Re] 


^  0.332 

Pr^Rel 


(5) 


where  Pr  is  Prandtl  Number  and  Re^  is  local  Reynolds  ITumber 
defined  by 


Re 


X 


/>e%^ 

Me 


(6) 


Note  that  yL4  is  fluid  viscosity  and  the  subscript  e  denotes 
boundary  layer  edge  conditions. 

For  real  gases  at  high  velocity  and  temperature,  the 
properties  do  vary  with  temperature.  They  may  be  assumed 
independent  of  pressure  as  long  as  excessively  high  pressure 
and  dissociation  are  avoided,  under  these  circumstances, 
and  as  long  as  c^  can  be  assumed  constant,  the  preceding 
relations  describing  heat  transfer  can  still  be  used  with 
good  accuracy  when  the  fluid  properties  are  calculated  at  a 
reference  temperature  T  ,  defined  as 

T*  =  Tg  +  0.5(T^  -  Tg)  +  U.22(Tj.  -  )  (7) 

by  Sckert  (Ref  2:10-13)« 

The  above  equations  may  be  used  to  solve  for  the  heat 
transfer  coefficient.  Prom  the  Nusselt  Number, 

h^  =  0.332  Pr%e^k/x  (8) 

Using  the  definition  of  Prandtl  iTuraber  as 


(9) 


the  heat  transfer  coefficient  can  be  reduced  to  the  expres¬ 


sion, 


h_  =  0.332 
^  Pr%  X 


Based  on  Eq.  (2),  definition  of  Stanton  llunber,  the  same  rela¬ 
tion  for  heat  transfer  coefficient  is  derived.  Ho^irever, 
other  definitions  of  Stanton  Number  are  given  in  various 
sources,  including 


3^x  = 


pc  U.0 


In  all  cases,  the  subscript  infinity,  oo,  refers  to  freestream 
conditions.  Based  on  these  varied  definitions  of  Stanton  Num¬ 
ber,  using  ideal  gas  calculations  of  fluid  properties,  several 
definitions  of  heat  transfer  coefficient  result,  prom  the 
Stanton  Number  of  Eq  (2), 

h^  =  0.332  ^  (10) 


Using  Eq  (11)  to  define  Stanton  Number, 

_  ,  .  ..2..n  i 


h^  =  0.332  _fp  T 

1 
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Defining  Stanton  Number  as  in  Eq  (12), 


(14) 


Eqs  (10),  (13) »  and  (14)  may  be  rewritten  for  the  high  speed 
case  by  substituting  for  local  density^  and  viscosityyt^with 
values  p  and/t  calculated  according  to  the  reference  tern- 
perature,  T  ,  defined  by  Eckert  as  in  Eq  (7).  The  resulting 
equations  are 


h^  =  0,332 

(15) 

Pr^j 

-  i 

h^  =  0,332 

c 

P. 

1 

(16) 

h^  =  0,332 

1 

r  P  2 

/fiw  r 

(17) 

X. 

Pr^l 

respectively.  In  addition,  Hayes  and  Neumann  (Ref  3:2)  in¬ 
troduce  the  empirical  calculation  of  the  ratio  of  heat  trans¬ 
fer  coefficients  as 


(18) 


The  above  equations  will  subsequently  be  used  to  evaluate 
reference  conditions. 


Method  of  Superposition 

The  method  of  wall  temperature  superposition  can  be  used 
to  develop  heat  transfer  solutions  for  boundary  layer  flow 
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over  a  plate  with  arbitrarily  varying  surface  temperature  due 
to  the  linearity  of  the  boundary  layer  energy  equation,  ac¬ 
cording  to  Kays  (Ref  4:218-219).  Assuming  constant  velocity 
flow  of  a  constant  property  fluid  over  a  flat  plate,  with 
laminar  flov;  only,  the  heat  transfer  coefficient  is  given  as 

h(§,x)  =  0.332  I  pr'^Rel  [l  -{^  (19) 


where  %  is  defined  as  the  unheated  plate  length.  Defining 
0(i»x,y)  as 

rp  q, 

0(f,x,y)  =  -2 -  (20) 

To  -  Too 

where  is  the  plate  surface  temperature.  Too  is  the  stream 
temperature  beyond  the  boundary  layer,  and  T  is  local  tem¬ 
perature,  a  solution  to  the  boundary  layer  energy  equation 
with  arbitrary  surface  temperature  distribution  is  given  by 

+  6  [l  -  e(1i»x,y)]  (21) 


The  heat  flux  from  the  wall  surface  is 


(22) 


Substituting  for  the  differential  with  Sq  (25)  yields 
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k 

(23) 

For  the  isothermal  wall  case,  with  wall  temperature  unequal 
to  freestream  temperature, 

4”  =  ke^($,x,0)  (T^  -  T«)  (24) 

By  convective  heat  transfer, 

4"  =  HHq  -  T^)  (25) 

Solving  for  the  heat  flux  in  terms  of  heat  transfer  coefficient 
and  conductivity  yields 

ey(t,x,0)  =  (26) 

Substituting  3q  (26)  into  Eq  (23)  gives  the  heat  flux  du-c  to 
an  arbitrary  wall  temperature  distribution  as 

where  the  heat  transfer  coefficient  is  defined  by  Eq  (23) 

(Ref  5:218-219)»  The  local  heat  transfer  coefficient  is  de¬ 
rived  directly  from  qJJ  as 

=  4”  (28) 

(Tq  -  To.) 

This  method  may  be  adapted  for  the  hypersonic  velocity  con¬ 
ditions  by  using  the  reference  temperature  defined  by  Eq  (7) 
to  calculate  the  fluid  properties,  and  by  substituting  the 
recovery  temperature  Tj.  for  the  freestream  temperature. 
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lov/  Temperature  Viscosity 

The  standard  relationship  used  to  calculate  viscosity 
is  that  given  by  Sutherland's  Law; 


_  2.27x10"^  T^ 
^  “  T  +  198.6 


lb-sec 

ft* 


(29) 


According  to  Piore  (Ref  5)»  this  may  not  be  valid  for  hyper¬ 
sonic  v.'ind  tunnel  applications.  In  hypersonic  tunnels,  the 
freestream  temperature  generally  falls  in  a  range  betv/een 
T  =  30  R  and  T  =  200  R  ,  depending  upon  the  ambient  tem¬ 
perature  of  the  gas  and  the  design  liach  mimber.  As  a  result, 
the  freestream  temperature  is  outside  the  demonstrated  range 
of  applicability  of  Sutherland's  Law,  as  an  alternative  cal¬ 
culation,  the  viscosity  can  be  defined  as 


2,32x10^  T^  lb-sec 

M.  =  —  , - r  — 

1  +/  220  ' 


(30) 


\Tx10 


This  relationship  provides  good  correlation  to  experimental 
data  for  lov/  temperature  flov/.  In  addition,  calculations  of 
viscosity  for  temperatures  above  T  =  300  R  by  this  low 
temperature  relationship  correspond  well  v/ith  viscosity  cal¬ 
culations  by  Sutherland's  Lav/  (Ref  6), 
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III.  Test  Data 


As  'indicated  in  the  introduction,  a  number  of  wind  tunnel 
tests  were  conducted  using  thin  skin  stainless  steel,  HRSI, 
and  PRSI  test  articles.  Prom  these,  three  runs  for  HRSI  and 
thin  skin  materials  were  chosen  to  yield  data  at  varied  de¬ 
flection  angles.  In  each  case,  data  vras  extracted  for  de¬ 
flection  angles  of  dX=  3,  10,  and  14  degrees  .  The  selec¬ 
tion  of  the  runs  used  for  each  case  v/as  based  on  injection  at 
high  deflection  angle  to  promote  rapid  temperature  response, 
and  on  length  of  run  at  the  desired  deflection  angle.  For 
HRSI  material,  experimental  runs  chosen  were  run  4  for 
0(  =  3  degrees  ,  run  20  for  o(  =  10  degrees  ,  and  min  19  for 
OC=  14  degrees  .  For  thin  skin  stainless  steel,  selected 
runs  were  run  58  for  oC  =  3  degrees  ,  run  57  for  oC=  10  de¬ 
grees  ,  and  run  46  for  <X=  14  degrees  •  In  addition,  data 
for  PRSI  material  runs  was  taken  from  runs  29  and  33  for 
=  3  degrees.  The  test  sequence  for  the  HRSI  cases  is  given 
in  Table  I,  for  the  PRSI  cases  in  Table  II,  and  for  the  thin 
skin  stainless  steel  cases  in  Table  III.  For  the  HRSI  runs, 
a  time  history  of  the  temperature  reading  on  thermocouple 
T/C  No.  5  is  given  in  pig  6,  with  location  given  in  Pig  1. 

For  the  PRSI  runs,  the  temperature  history  is  given  in  Pig  7 
for  thermocouple  t/C  No.  4,  with  location  given  in  pig  2.  por 
the  thin  skin  runs,  the  temperature  history  for  thermocouple 
T/C  no.  4  is  given  in'Tj’ig  8,  with  location  given  in  Pig  3. 

The  thermocouples  were  chosen  to  give  temperature  history 
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data  at  the  material  surface.  Pig  6  reveals  that  the  HRSI 
test  article  levels  out  at  a  temperature  of  T  =  705  R  for 
a  deflection  angle  of  0C=  3  degrees  before  the  angle  is 
raised  to  oC=  14  degrees  ,  in  run  4»  Similar  indications 
are  gained  from  run  20,  where  T  =  760  R  for  oC=  10  degrees 
and  from  run  19,  where  T  =  810  R  for  cL-  14  degrees  • 

Prom  these,  it  would  appear  that  thermal  equilibrium  is  a- 
chieved  in  each  case.  This  equilibrium  point  gives  the  maixi- 
mum  wall  temperature  achieved  at  the  thermocouple  location 
used.  A  similar  equilibrium  point  is  seen  in  the  time  history 
data  for  PR3I  material  in  run  33,  Pig  7.  This  case  shov;s 
thermal  equilibrium  for  PP3I  material  at  a  deflection  angle 
of  0C=  3  degrees  to  be  T  =  700  R  .  Por  run  29,  it  appears 
that  equilibrium  is  approached  more  gradually  and  is  not  a- 
chieved  before  the  test  run  is  ended.  As  a  result,  PRSI  sur¬ 
face  heating  data  were  taken  from  run  33  for  comparison. 

At  the  thermal  equilibrium  points  for  the  HHSI  and  ?RSI 
materials,  data  readings  were  extracted  for  all  thermocouple 
locations.  It  was  found  that  the  HH3I  thermocouple  t/c  Ho,  6 
was  not  functioning  during  any  of  the  runs  made,  as  it  main¬ 
tained  a  constant  temperature  reading.  The  location  of  this 
thermocouple  is  given  in  Pig  1,  As  a  result,  only  one  surface 
data  point  was  available  for  HR 31  heating  results,  Por  PRSI 
material,  run  33  was  used  to  extract  a  surface  temperature 
distribution  at  a  deflection  angle  of  oC=  3  degrees  .  The 
temperature  readings  for  each  location  are  given  in  pig  9. 
prom  this  data,  the  P7.SI  surface  tem::)erature  is  revealed  to 
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vary  by  only  a  fev/  degrees  above  and  below  a  constant  of 
T  =  690  R  •  This  is  a  strong  indication  of  constant  tem¬ 
perature  over  the  insulating  material.  Since  the  thermal 
design  properties  of  the  HRSI  material  are  similar  to  those 
of  the  pnsi  material,  the  surface  temperature  of  the  HRSI 
material  vrill  be  assumed  constant  streamwise,  at  the  value 
given  by  HRSI  thermocouple  t/C  Ko.  5  for  each  deflection 
angle. 

The  thermocouple  data  for  the  thin  skin  temperature  dis¬ 
tribution  at  each  deflection  angle  were  taken  from  runs  58, 

57,  and  46,  as  indicated  e^'nlier,  A  temperature  time  history 
for  the  thin  skin  article  is  given  in  Pig  8  for  these  runs. 

It  is  seen  here  that  the  surface  temperature  increases 
throughout  the  test  runs  for  the  thin  skin  cases,  since  no 
equilibrium  was  achieved,  the  surface  temperature  distribution 
data  v;ere  arbitrarily  taken  at  the  point  of  maximum  value, 
near  the  end  of  each  run.  The  distribution  for  deflections 
of  Oi=  10,  and  14  degrees  are  given  in  Pigs  9,  10,  and 
11,  respectively.  For  each  case,  the  thin  skin  stainless 
steel  surface  temperature  gradually  decreases  v^ith  streamwise 
location.  This  decline  is  a  very  small  percentage  of  the 
surface  temperature  itself,  so  the  thin  skin  surface  is  as¬ 
sumed  isothermal  at  a  median  wall  temperature.  Since  no  data 
was  available  for  the  wedge  surface  temperature  ahead  of  the 
test  article  location,  the  initial  wall  temperature  was  as¬ 
signed  the  thin  skin  value  for  the  given  deflection  angle, 
and  was  assumed  constant. 
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It  should  be  noted  that  several  underlying  assumptions 
made  here  did  not  allov;  duplication  of  the  test  conditions  in 
the  numerical  calculation.  The  HRSI  thermocouple  actually 
lies  beneath  the  radiative  coating  on  the  surface,  and  may 
not  reflect  the  actual  surface  temperature.  The  coating  sur¬ 
face  temperature  may  be  somewhat  higher.  The  higher  deflec¬ 
tion  angle  readings  for  the  HRSI  may  not  be  at  equilibrium 
but  may  be  only  approaching  that  condition.  A  higher  tempera¬ 
ture  at  equilibrium  would  then  be  expected.  Since  the  wedge 
surface  temperature  is  unknown,  it  may  not  be  at  the  tempera¬ 
ture  of  the  thin  skin  plate,  indeed,  the  wedge  is  a  much 
larger  heat  sink  than  the  thin  skin  plate  in  terms  of  volume. 
As  such,  its  surface  temi)erature  should  not  rise  as  quickly 
as  the  thin  skin  surface  temperature,  but  should  probably 
remain  close  to  the  initial  ambient  temperature  of  the  wedge. 

A  possible  indication  of  this  temperature  difference  betv;een 
tiie  thin  skin  and  wedge  stainless  steel  surfaces  is  seen  in 
the  heat  transfer  coefficient  data  given  by  Hayes  (Ref  8)  and 
displayed  in  Fig  12.  Note  that  this  data  reveals  that  the 
heat  transfer  coefficient  actually  rises  to  a  maximum  at  the 
location  x  =  9  in  ,  and  then  declines.  For  an  isothermal 
v;all,  the  distribution  of  local  heat  transfer  coefficient 
should  decrease  gradually  along  this  entire  length,  as  indi¬ 
cated  in  Fig  14,  v/ith  no  rise  like  th.at  seen  in  the  test  data. 
As  a  result,  some  v;all  temperature  discontinuity  may  be  pres¬ 
ent,  Note,  hov.'ever,  that  this  rise  in  the  heat  transfer  co¬ 
efficient  may  also  be  caused  by  bov/ing  of  the  heated  thin  skin 
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plate  as  a  result  of  expansion.  Higher  heating  would  he  seen 
on  the  ^forward  part  of  the  plate,  with  somewhat  lower  heating 
on  the  “aft  section  for  this  case  as  v/ell.  Also,  constant 
temperature  on  the  viedge  ahead  of  the  test  article  location 
was  assumed,  but  could  not  be  verified.  Finally,  the  thin 
skin  readings  were  taken  at  their  maximum  values,  v;hich  may 
not  correspond  to  the  same  time  coordinate  as  the  HR3I  equili¬ 
brium  readings.  The  indication  of  these  differences  is  that 
the  initial  vrall  temperature  should  probably  be  lov:er  than 
the  thin  skin  values  assumed,  the  final  wall  temperature  may 
be  somewhat  higher,  and  the  step  size  is  most  likely  larger 
than  that  taken  from  the  test  data. 
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IV.  Reference  Conditions 


Reference  calculation 

In  studying  the  variation  of  heating  on  the  HR3I,  PRSI, 
and  thin  skin  test  articles  due  to  deflection  angle  and  wall 
temperature  variation,  boundary  layer  characteristics  become 
very  important.  The  heat  transfer  to  the  surface  depends  on 
the  thermal  transport  properties  of  the  boundary  layer.  The 
convective  heat  flux  is  defined  as  the  proiuct  of  the  tempera¬ 
ture  potential  and  the  heat  transfer  coefficient.  The  tem¬ 
perature  potential  is  the  difference  between  the  wall  tempera¬ 
ture  and  the  freestream  temperature  for  low  speed  flow,  or 
between  the  wall  temperature  and  recovery  temperature  for  high 
speed  flov^.  As  such,  the  temperature  potential  is  independent 
of  boundary  layer  properties.  As  a  result,  boundary  layer 
heat  transport  dependence  on  v/all  temperature  variation  must 
be  reflected  in  the  heat  transfer  coefficient  distribution. 

In  ord^r  to  enh-nce  the  comparison  of  results,  and  elimi¬ 
nate  location  dependence  of  the  heat  transfer  coefficient  in 
the  comparison,  a  reference  distribution  at  zero  deflection 
was  required.  The  reference  value  at  location  x/l  =  0,575 
is  of  particular  significance  since  HRSI  data  are  given  for 
this  location.  This  value  was  calculated  by  the  varied  defini¬ 
tions  of  heat  transfer  coefficient  according  to  Eckert  flat 
plate  and  high  speed  flat  plate  theory.  At  a  deflection  angle 
of  zero  degrees,  the  boundary  layer  edge  values  and  freestream 
values  of  density,  velocity,  and  viscosity  are  identical,  and 


Sqs  (10),  (13),  and  (14)  yield  the  same  results.  The  reference 
heat  transfer  coefficient  calculated  "by  Eckert  flat  plate 
theory  Is  =  1*12  x10“^  Btu/ft‘^s-R  ,  based  on  the  v^ind 

tunnel  freestream  operating  conditions,  Prandtl  iTumber  is 
assumed  constant  at  pr  =  0,72  ,  and  the  specific  heat  is 
assumed  constant  at  c^  =  6006  ft^/s^  .  These  values  v;ill 
be  assumed  constant  in  all  subsequent  calculations.  According 
to  Eckert  high  speed  flat  plate  theory,  the  fluid  properties 
are  determined  according  to  the  reference  temperature  defined 
in  Eq  (7).  At  zero  degree  deflection,  Eqs  (15),  (16),  and 
(17)  yield  the  same  result  for  heat  transfer  coefficient, 
h^ef  =  xIO""^  3tu/ft^s-R  . 

As  indic'^.ted  in  the  Theory  section,  Sutherland's  Viscosity 
Lav;  may  not  bo  valid  at  lov;  freestream  temperatures.  An  al¬ 
ternate  form  given  by  Eq  (30)  can  be  used  to  estimate  a  refer¬ 
ence  heat  transfer  coefficient.  The  variation  between  this 
relationship  and  that  given  by  Sutherland's  Lav;  is  negligible 
above  a  temperature  of  T  =  500  R  ,  as  seen  in  Pig  13,  Be¬ 
low  T  =  300  R  ,  there  is  significant  divergence,  v^ith  the 
low  temperature  calculation  yielding  viscosity  values  greater 
than  those  given  by  Sutherland's  Lav;  at  the  same  temperature, 

A3  a  result,  v;hen  applied  to  Eq  (10),  the  lov;  temperature  cal- 

—4 

culation  reveals  a  ncv;  reference  value  of  h  -  =  1,29  x10 

ref 

2 

Btu/ft  s-R  .  pdien  applied  to  Eq  (15),  hov;ever,  the  reference 
calculation  is  unchanged  because  the  fluid  properties  are 
determined  at  the  reference  temperature,  T*  =  670  R  .  As  a 
result,  the  lov;  temperature  viscosity  problem  is  only  consi- 
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derable  when  using  the  basic  flat  plate  theory. 

Thel  boundary  layer  solution  by  numerical  technique  is  ac- 
compli^ed  by  the  boundary  layer  program  referenced  in  Appen¬ 
dix  B*  Reading  the  freestream  conditions  of  Mach  Number 
temperature  T«a,  Prandtl  Humber  Pr  (as  assumed  earlier),  and 
specific  heat  ratio  (assumed  as  r=  1.4-  in  all  cases), 


the  program  yields  the  local  Stanton  Number  3t^  from  which  the 
local  heat  transfer  coefficient  is  derived.  The  reference 
value  given  by  this  method  is  =  9»5  x10“^  Btu/ft^s-R 


for  a  wall  temperature  of  =  530  R  ,  and 
x10~^  Btu/ft^s-R  for  a  wall  temperature  of 


\ef  =  9.55 


=  560  R 

w 


A  reviev;  of  the  program  calculation  revealed  that  the  recovery 


factor  v/as  not  accounted  for  in  estimating  the  temperature 


potential,  in  hypersonic  calculations,  this  recovery  factor 


r  is  generally  assumed  as  r  =  0.9  .  According  to  Holman 
(Ref  8:213)  as  vrell  as  others,  for  laminar  flow,  r  =  J"pr  , 


Assuming  Pr  =  0.72  ,  then  r  =  0.85  .  Recalculation  of  the 


program  value  of  reference  heat  transfer  coefficient  by  each 
of  these  yields  =  1.21  x10“'^  Btu/ft^s-R  for  r  =  0.85  , 

and  =  1.11  xIO”^  Btu/ft^s-R  for  r  =  0.9  •  The  latter 

of  these  is  consistent  with  the  theoretical  values.  All  re¬ 


maining  estimations  'ill  be  based  on  a  reference  of  = 

1.1  x10“^  Btu/ft'^s-R  as  it  is  consistent  with  theory  and  with 
the  original  test  reference.  Additionally,  a  recovery  factor 
of  r  =  0.9  will  be  used  in  all  subsequent  calculations, 
both  numerical  and  analytical.  This  is  consistent  with  the 
original  data  reduction  also. 
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An  important  consideration  of  the  varying  surface  tem¬ 
perature,  or  nonisothermal  wall  effect,  is  the  streamwise 
variation  of  heat  transfer  coefficient  beyond  the  temperature 
step.  The  distance  downstream  over  vrhich  this  effect  is  felt 
must  be  v;ell  defined  in  order  to  accurately  evaluate  heat 
transfer  beyond  the  wall  temperature  discontinuity.  The 
variation  of  the  heat  transfer  coefficient  for  an  isothermal 
wall  is  then  needed  for  comparison.  According  to  Eckert  high 
speed  theory,  Eq  (15),  this  variation  is  shovm  in  Pig  14. 

As  indicated,  the  heat  transfer  coefficient  decreases  hyper- 
bolically  from  the  leading  edge  and  levels  out  downstream, 

A  comparison  is  given  in  Pig  15  betvreen  heat  transfer  coeffi¬ 
cient  calculations  according  to  Eckert  high  speed  theory  and 
the  boundary  layer  program  for  both  Sutherland  and  lov:  tem¬ 
perature  viscosity.  As  shovm,  the  variation  of  heat  transfer 
coefficient  is  consistent  v;ith  that  indicated  in  Pig  14. 

Note  that  these  values  correspond  to  the  isothermal  vmll  case 
at  zero  deflection. 

Isothermal  'jail  Calculations 

In  the  reference  calculations  above,  the  various  flat 
plate  theory  relations  collapsed  to  one  expression,  as  did 
the  high  speed  theory  equations,  since  the  boundary  layer  edge 
values  match  the  freestream  values  v/ith  no  deflection  of  the 
plate.  At  a  given  deflection,  however,  an  oblique  shock  occurs 
at  the  leading  edge  of  the  plate  and  the  boundary  layer  edge 
conditions  are  no  longer  the  freestream  properties.  The  proper¬ 
ties  at  the  boundary  layer  edge  are  derived  from  the  oblique 
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shock  relations  "based  on  freestream  conditions.  Using  these 
values,-  the  various  relations-derived  for  flat  plate  theory, 

Sqs  (10),  (15),  and  (14),  are  revealed  to  diverge  drastically. 
Similarly,  the  results  gained  fron  the  definitions  of  Sqs  (15)- 
(17)  for  high  speed  flow  bear  some  disagreement.  Using  the 
boundary  layer  edge  conditions,  the  boundary  layer  program 
revealed  isothermal  v;all  results  which  proved  to  be  very  con¬ 
sistent  with  the  Eckert  high  speed  calculation  of  Eq  (15). 

Due  to  this  correlation  betv/een  theory  and  the  numerical  solu¬ 
tion,  Eq  (15)  has  been  found  to  be  the  correct  expression  of 
the  Eckert  theory,  v;ith  the  definition  of  Stanton  Number  given 
in  Eq  (2).  The  corresponding  low  speed  expression  is  that 
given  by  Eq  (10).  The  results  of  Bqs  (15)  and  (14)  for  low 
speed  flow,  and  of  Eqs  (16)  and  (17)  for  high  speed  flow  were 
found  to  be  unacceptable  when  compared  to  the  numerical  solu¬ 
tion.  The  correspondence  between  theory  and  numerical  solu¬ 
tion  is  given  in  ?ig  l6.  Also  of  note  are  the  exactness  of 
the  approximation  given  by  Eq  (18)  and  the  divergence  betvreen 
the  low  speed  and  high  speed  theoretical  solutions  at  increas¬ 
ing  deflection  angles.  Again,  the  numerical  solution  verifies 
the  high  speed  theory  in  application  to  this  problem.  The 
correspondence  between  theory  and  numerical  analysis  lends 
further  support  to  the  correction  of  the  program  for  a  recovery 
factor  of  r  =  0.9  . 

The  boundary  layer  program  was  modified  as  indicated  in 
Appendix  B  to  solve  the  exact  equations  using  the  lov;  tempera¬ 
ture  viscosity  calculation  of  Eq  (30).  Using  the  boundary 
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layer  edge  condiuions,  the  program  calculations  revealed  that 
the  heart  transfer  coefficient  for  the  isothermal  v/all  and  low 
temperature  viscosity  did  not  diverge  significantly  from  the 
Sutherland's  Law  calculations  hut  agreed  very  closely.  Cal¬ 
culation  vith  So  (15),  Eckert  high  speed  theory,  agreed  pre¬ 
cisely  with  the  boundary  layer  program  for  this  case  also, 
pig  17  reveals  the  close  correlation  between  theory  and  nu¬ 
merical  solution  once  again.  Note  the  divergence  between 
high  speed  and  low  speed  theory  for  increasing  angle  of  de¬ 
flection,  A  comparison  of  program  data  betv/een  Pigs  16  and 
17  reveals  the  close  correlation  for  the  tv;o  different  vis¬ 
cosity  calculations.  This  further  supports  the  high  speed 
theory  based  on  Eckert's  reference  temperature.  The  high  speed 
theory  and  the  boundary  layer  isothermal  wall  calculations 
v;ill  be  the  basis  of  comparison  for  the  nonisothermal  v;all 
boundary  layer  calculations. 


24 


V*  Numerical  Results 

Nonisothermal  Wall 

Using  the  boundary  layer  program  modified  for  a  wall  tem¬ 
perature  step,  as  indicated  in  Appendix  B,  the  numerical  sol¬ 
ution  of  the  boundary  layer  equations  can  be  derived  for  a  con¬ 
stant  temperature  v;all  with  a  step  change.  This  condition  is 
the  same  as  that  assumed  from  the  v/ind  tunnel  test  data.  The 
initial  temperature  of  the  wall  is  that  gained  from  the  thin 
skin  test  data  at  each  deflection  angle.  For  a  deflection 
angle  of  <\=  3  degrees  ,  a  thin  skin  wall  temperature  of 
T^  =  600  R  was  extracted,  similarly,  for  oC=  10  degrees  , 

Tw  =  ^20  R  ,  and  for  =  14  degrees  ,  T^  =  640  R  .  These 
readings  are  also  the  basis  of  the  isothermal  wall  calcula¬ 
tions  in  the  boundary  layer  program,  as  revealed  in  the  pre¬ 
vious  section.  As  given  earlier,  the  wall  temperature  over 
the  HR3I  material  ^ assumed  constant)  is  the  new  wall  tempera¬ 
ture  after  the  step  increase.  For  the  deflection  oC=  3  de¬ 
grees  ,  the  wall  temperature  after  the  step  is  T^^  "  ^  • 

Similarly,  for  0(=  10  degrees  ,  ^  OC  = 

14  degrees  ,  T^2  =  810  R  •  Again,  the  ratio  of  heat  trans¬ 

fer  coefficient  to  the  ref'rence  value  is  calculated  for  a 
strearawise  location  of  x/l  =  0,575  , 

Based  on  the  given  initial  and  final  wall  temperatures, 
the  derived  values  of  the  ratio  h/h^^^  for  the  step  change 
in  wall  temperature  are  given  in  Fig  18,  The  isothermal  case, 
based  on  thin  skin  v^all  temperature  at  each  angle  of  deflec- 
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tion,  is  given  for  comparison.  Note  that  the  nonisothermal 
wall  result  shov/s  an  increasing  divergence  from  the  isothermal 
case  with  increasing  deflection  angle.  At  a  deflection  of 
OC  =  3  degrees  ,  the  nonisothermal  case  diverges  from  the  iso¬ 
thermal  by  twelve  percent  v/hile  the  divergence  is  nearly 
tv;enty  percent  at  an  angle  of  0C=  14  degrees  .  Also  given 
for  comparison  are  the  ^lEATUST  data  results  based  on  HR3I 
temperature  data,  and  an  equilibrium  calculation  based  on 
HR'^I  temperature  at  a  deflection  angle  of  0C=  3  degrees  . 

In  this  figure  it  is  seen  that  the  nonisothermal  effects  ac¬ 
count  for  only  about  one-fourth  of  the  discrepancy  between  the 
HEAT^^ST  data  and  the  isothermal  v/all  case  at  a  deflection 
angle  of  OC=  3  degrees  .  However,  at  (X=  14  degrees  , 
the  HEATEIT  result  is  identical  to  the  nonisothermal  case. 

For  each  angle  of  deflection  studied,  the  variation  of 
heat  transfer  coefficient  v;ith  streanv;ise  location  is  derived. 

As  seen  in  Pigs  19,  20,  and  21,  for  oC=  3,  10,  and  14  de¬ 
grees  ,  respectively,  the  isothermal  data  shov/s  the  same  trends 
as  indicated  by  the  reference  case  in  pig  15.  Por  each  de¬ 
flection  angle,  the  streamv/ise  location  of  the  wall  tempera¬ 
ture  step  is  x/l  =  0.467L  v/here  L  is  the  plate  length  and 
is  given  to  be  L  =  1.25  ft  .  According  to  the  boundary 
layer  program,  the  step  increase  in  temperature  creates  a 
discontinuity  in  the  heat  transfer  coefficient,  after  which 
it  rises  sharply,  levels  out  quickly,  and  gradually  decreases 
tov/ard  the  trailing  edge  of  the  plate  ^  as  sho\'m.  After  the 
discontinuity  then,  the  heat  transfer  coefficient  approaches 
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the  condition  of  the  isothermal  case,  but  sees  a  discrete 
decrease  in  magnitude.  The  same  is  seen  for  all  three  de¬ 
flection  angles  discussed.  Note  that  the  drop  in  heat  trans¬ 
fer  coefficient  is  sharper  v/ith  increasing  angle  of  deflection, 
and  recovery  toward  the  isothermal  case  is  slower.  Lastly, 
the  magnitude  of  the  difference  between  isothermal  and  non- 
isothermal  cases  is  increasing  with  greater  deflection  angle, 
from  about  six  percent  of  the  isothermal  value  at  the  trailing 
edge  at  OC=  3  degrees  to  nine  percent  at  oC=  14  degrees  • 

As  indicated  earlier,  the  change  in  results  due  to  the 
low  temperature  viscosity  relationship  is  very  small  since 
the  case  considered  is  at  very  high  speed  and  thus  has  a  high 
reference  temperature.  This  is  verified  once  again  in  calcu¬ 
lating  the  nonisothermal  v;a.ll  case.  Considering  the  values 
shovm  in  pig  22,  there  is  very  little  difference  betv;een  the 
isothermal  values  here  and  those  shown  in  pig  21.  The  same 
is  true  for  the  nonisothermal  v/all  values  for  the  ratio  of 
of  heat  transfer  coefficients.  The  same  trends  are  sho\im  as 
angle  of  deflection  increases,  v^ith  the  divergence  grovring 
betv;een  isothermal  and  nonisothermal  cases.  Consequently, 
it  has  been  shov;n  through  graphical  results  of  the  isothermal 
and  nonisothermal  wall  calculations,  a.nd  through  the  variation 
of  heat  transfer  coefficient  with  location,  that  the  lov;  tem¬ 
perature  viscosity  calculation  has  shovm  no  difference  in  the 
trends  of  behavior  v;ith  location,  deflection  angle,  or  step 
presence,  and  very  little  difference  in  the  magnitude  of 
heat  transfer  coefficient  ratio  in  any  of  these  cases. 
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step  Size  Variation 

In  reducing  the  test  data,  several  considerations  indi¬ 
cate  that  the  step  size  and  wall  temperatures  assumed  may  not 
he  accurate.  The  effect  of  a  lov;er  wedge  temperature,  coupled 
with  a  higher  IIHSI  surface  equilibrium  temperature,  would 
cause  a  significsjit  increase  in  step  size,  with  associated' 
boundary  layer  effects.  A  larger  step  size  should  certainly 
have  an  effect  on  the  magnitude  and  recovery  distance  of  the 
heat  transfer  coefficient.  Thus,  understanding  the  effects 
of  a  larger  step  would  be  beneficial  in  assessing  the  validity 
of  conclusions  drawn  from  other  results.  In  addition,  in¬ 
creasing  the  step  size  makes  the  test  results  more  versatile, 
and  more  broadly  applicable  to  a  range  of  considerations. 

To  investigate  the  effect  of  increasing  the  step  size, 
the  step  used  previously  in  each  deflection  angle  case  was 
incremented.  The  heat  transfer  coefficient  variation  was  de¬ 
rived  for  a  step  of  twice  the  original  step,  four  times  the 
original  step,  and  tv^ice  the  initial  wall  temperature.  This 
v;as  accomplished  for  each  deflection  angle.  The  results  for 
a  deflection  angle  of  0(  =  3  degrees  are  plotted  in  Pig  25, 
for  oC=  10  degrees  in  Pig  24,  and  for  oC  =  14  degrees  in 
Pig  25.  Por  the  case  of  deflection  angle  0C=  5  degrees  , 
the  heat  transfer  coefficient  sees  a  sharp  discontinuity  for 
all  cases.  This  discontinuity  grows  with  step  size.  Por 
smaller  steps,  the  heat  transfer  coefficient  rises  sharply, 
peaks,  and  decreases  gradually  toward  the  trailing  edge. 

Por  increasing  step  size,  the  minimum  value  is  no  longer  seen 
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immediately,  but  comes  more  gradually  as  heat  transfer  coef¬ 
ficient  falls  more  slowly,  then  recovers  more  slowly.  For 
the  step  size  of  four  times  the  original  step  size,  a  peak  is 
not  reached  until  x/l  =  0,8  ,  For  still  larger  steps,  the 

peak  may  not  be  reached.  Generally,  reaction  of  heat  trans¬ 
fer  coefficient  is  slov/er  for  increasing  step  size,  while 
minimum  value  is  smaller  and  magnitude  dovmstream  is  smaller. 
This  is  also  seen  to  be  true  for  the  same  increments  of  step 
size  at  deflections  of  a  =  10  and  14  degrees  ,  In  addition, 
note  that  reaction  and  recovery  are  slov/ing,  and  magnitude 
decreasing  as  the  deflection  angle  increases,  as  noted  earlier 
and  emphasised  v.'ith  this  data. 

The  size  of  the  step  can  also  be  changed  by  lov;ering  the 
initial  v/all  temperature.  This  is  consistent  v;ith  the  lo\reT 
temperature  expected  on  the  stainless  steel  v;edge,  as  men¬ 
tioned  previously.  Indeed,  from  experience,  the  wedge  tem¬ 
perature  may  remain  very  close  to  its  initial  temperature, 

'H’or  the  deflection  angle  OC  =  3  degrees  ,  the  surface  tem¬ 
perature  initially  is  T^^  =  540  R  ,  and  at  o(=  14  degrees  , 
the  surface  temperature  is  T^^  =  550  R  ,  Using  these  ini¬ 
tial  wall  temperatures,  v/ith  the  HTiGI  final  temperature  for 
the  appropriate  deflection  angle,  yields  the  streamwise  dis¬ 
tribution  of  heat  transfer  coefficient  for  oC  =  3  and  14  de¬ 
grees  in  Figs  26  and  27,  re  selectively.  The  heat  transfer 
coefficient  ratio  for  each  deflection  angle  is  shown  in  Pig 
18,  These  values  reveal  that  at  a  deflection  angle  of 
a  =  3  degrees  ,  the  nonisothermal  wall  accounts  for  about 
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forty  percent  of  the  discrepancy,  while  at  a  deflection  of 
(X  =  14  degrees  ,  about  sixty  percent  is  accounted  for,  as 
seen  in  Pig  18.  These  results  require  some  qualification. 

The  initial  wall  temperature  used  here  is  the  lov/est  value 
possible,  before  the  wedge  is  subjected  to  any  heat  transfer. 
This  should  be  reasonable  for  the  lov;  angle  case,  but  is 
questionable  in  the  higher  angle  case,  v.'here  heating  is  sig¬ 
nificantly  higher.  Therefore,  the  reliability  of  the  lower 
angle  case  is  probably  better.  Comparing  these  values  of  the 
heat  transfer  coefficient  ratio  v/ith  pig  28,  the  wall  tem¬ 
perature  step  may  have  to  be  double  that  assumed  here  to  ac¬ 
count  for  the  entire  discrepr’jicy,  ’ihile  the  PPSI  surface 
temperature  may  be  higher  than  that  indicated  by  the  test 
data,  it  should  be  doubtful  that  it  will  be  tv:ice  that  in¬ 
dicated. 

The  variation  of  the  ratio  of  heat  transfer  coefficients 
with  increasing  step  size  and  angle  of  deflection  is  revealed 
in  pig  28,  The  isothermal  wall  case  is  given  for  comparison 
and  all  values  are  shown  for  a  location  of  x/l  =  0,575  • 

As  shown,  there  is  no  real  pattern  v/ith  increasing  step  size 
except  that  the  ma.gnitude  is  decreasing.  Note,  however,  that 
the  effect  of  recovery  location  is  important.  At  deflection 
of  OC  =  5  degrees  ,  all  but  the  largest  step  are  fairly  well 
recovered  by  x/l  =  0.575  .  Only  the  magnitude  of  the  largest 

step  is  very  small.  At  an  angle  of  0C=  10  degrees  ,  the 
step  of  four  times  the  original  step  is  at  its  minimum  value 
at  this  X  location,  while  the  largest  step  is  still  approach- 
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ing  minimum.  The  other  smaller  steps  are  v/ell  recovered. 
Finally,  at  CC=  14  decrees,  neither  of  the  tv^o  larger  steps 
has  reached  a  minimum,  while  the  two  smaller  steps  are  re¬ 
covering  at  x/l  =  0,575  .  This  serves  to  emphasize  the  de¬ 

lay  in  recovery  for  increasing  step  size  and  increasing  de¬ 
flection  angle,  as  v/ell  as  the  decrease  in  heat  transfer  co¬ 
efficient  for  increasing  step  size. 

Step  Lopg^-tion  Variation 

AS  an  additional  consideration,  the  effect  of  varying  the 
step  location  has  been  investigated.  The  results  are  shown 
in  pig  29,  This  set  of  calculations  was  conducted  using  a 
constant  temperature  step  /ith  a  constant  initial  wall  tem¬ 
perature  and  fixed  deflection  angle.  The  conditions  chosen 
were  at  a  deflection  angle  of  0C=  5  degrees  ,  v;ith  the  ini¬ 
tial  v/all  temperature  of  T^^  =  600  R  ,  and  a  final  wall 
temperature  of  T^2  =  ^  •  S'tep  locations  were  taken  as 

x/l  =  0,1,  0,25,  0,47,  0.6,  and  0,6  •  For  the  thermocouple 

location  x/l  =  0,575  ,  only  those  step  locations  upstream 

would  affect  the  heat  transfer  coefficient  seen.  However, 
from  the  results  it  is  seen  that  recovery  is  somev^hat  slower 
as  the  step  is  moved  downstream.  Additionally,  the  drop  in 
heat  transfer  coefficient  is  larger  ’..’ith  increasing  x  loca¬ 
tion,  Finally,  the  magnitude  of  the  deviation  from  the  iso¬ 
thermal  case  seems  to  increase  dov.’nstream,  although  it  doesn't 
become  very  large  with  the  given  step  size.  The  overall  ef¬ 
fect  is  to  incre.ase  recoverj'  distance  and  deviation  with  in¬ 
creasing  dov/nstream  step  location. 
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In  an  effort  to  verify  the  boundary  layer  program  results 
analytically  and  at  the  same  time  identify  the  trends  in  var¬ 
iation  of  heat  transfer  coefficient  v;ith  arbitrary  wall  tem¬ 
perature  variations,  a  secondary  means  of  estimating  the  non- 
isothermal  effect  v/as  desired.  An  analytical  method  of  esti¬ 
mating  the  heat  transfer  coefficient  v/ith  the  nonisothermal 
wall  v;hich  was  presented  in  the  section  on  Theory  is  the  super¬ 
position  technique.  This  method  was  applied  to  the  cases  of 
deflection  angle  OC  =  3,  10,  and  14  degrees  ,  with  the  cor¬ 
responding  initial  and  final  wall  temperatures.  As  used  in 
hypersonic  theory  calculations,  the  fluid  properties  v^ere  de¬ 
rived  from  the  reference  temperature  for  each  case,  and  the 
recovery  temperature  was  used  in  the  temperature  potential 
in  place  of  the  free stream  temperature. 

The  results  of  this  estimation  method  are  given  in  pigs 
50,  51,  and  52,  for  the  deflection  angles  of  0C=  3»  10,  and 
14  degrees  ,  respectively.  The  data  suggests  several  con¬ 
clusions,  The  change  in  heat  transfer  coefficient  from  the 
isothermal  distribution  increases  with  deflection  angle. 

Also,  the  superposition  estimation  does  demonstrate  great  dis¬ 
continuity  effects.  4  sharp  drop  and  recovery  are  shov/n  in 
each  case,  increasing  in  magnitude  with  deflection  angle. 

This  effect  is  similar  to  that  demonstrated  by  the  numerical 
solution. 

The  results  yielded  by  the  numerical  and  analytical 


methods  are  compared  in  pig  33.  The  methods  of  numerical 
solution  and  superposition  are  compared  to  the  isothermal  wall 
case  calculated  by  the  boundary  layer  program  and  Eckert  high 
speed  theory.  The  values  shown  represent  the  ratio  of  heat 
transfer  coefficient  at  the  streamwise  location  z/l  =  0,575  » 

for  varying  deflection  angle.  As  seen,  the  numerical  solu¬ 
tion  yielded  the  greatest  divergence  from  the  isothermal  wall 
calculations.  This  divergence  ranges  from  a  value  of  12  to 
19  percent  of  the  isothermal  value.  By  comparison,  the  range 
of  divergence  shovm  for  the  superposition  method  is  between 
6  and  13  percent.  The  result  is  that  the  analytical  results 
given  by  wall  temperature  superposition  reveal  consistent 
trends  in  the  heat  transfer  coefficient  change  due  to  the  wall 
temperature  discontinuity,  at  varied  deflection  angles.  How¬ 
ever,  the  magnitude  of  this  effect  is  consistently  smaller  ac¬ 
cording  to  superposition  t-.-.n  that  indicated  by  numerical 
solution.  The  values  derived  from  each  method  do  form  a  small 
range  of  values  in  v/hich  the  nonisothermal  v;all  effect  may 
be  expected  to  fall,  v/ith  the  limits  of  this  range  formed  by 
the  superposition  estimation  at  the  upper  limit  and  the  nu¬ 
merical  solution  at  the  lovrer  limit. 


VII. 


Conclusions  and  Recommendations 


Conclusions 

Experimental  Data.  The  test  data  revealed  several  uncer¬ 
tainties  that  left  the  actual  step  size  unknovm.  At  a  deflec¬ 
tion  angle  of  OC  =  3  degrees  ,  the  HRSI  test  article  reached 
equilibrium.  At  oC=  10  and  14  degrees  ,  this  condition  was 
not  achieved.  The  HRSI  temperature  used  for  these  cases  was 
lov^er  than  that  required  for  equilibrium.  The  thermocouple 
itself  was  not  at  the  surface  but  beneath  the  radiative  coat¬ 
ing,  The  associated  temperature  reading  may  have  been  lov/er 
than  the  actual  surface  temperature.  The  constant  surface 
temperature  assumption  for  the  HRSI  article  appears  valid. 

The  stainless  steel  v:edge  surface  temperature  was  not  neces¬ 
sarily  equal  to  the  appropriate  thin  skin  plate  surface  tem¬ 
perature,  This  is  strongly  indicated  by  the  heat  sink  argu¬ 
ment  and  by  the  thin  skin  heat  transfer  coefficient  data  cited. 
Constant  temperature  on  the  wedge  vras  assumed  but  could  not 
be  verified.  The  conclusions  from  this  data  is  that  initial 
wall  temperature  should  be  lower  than  the  thin  skin  data  used, 
HRSI  equilibrium  surface  temperature  should  be  higher,  espe¬ 
cially  for  the  higher  deflection  angle  cases.  The  overall 
step  was  significantly  larger  than  that  assumed  from  thin 
skin  and  HRSI  test  data. 

Reference  Conditions,  In  calculating  the  reference  value 
of  the  heat  transfer  coefficient,  Eckert  high  speed  flat  plate 
theory  proved  to  be  consistent  v/ith  the  boundary  layer  pro- 
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gram#  svaluation  of  program  results  considering  recovery 
temperature  showed  precise  correlation  between  numerical  and 
analytical  reference  solutions,  assuming  a  recovery  factor  of 
r  =  0.9  •  This  was  consistent  with  the  original  data  reduc¬ 

tion  as  well.  Low  speed  theory  diverged  drastically  at  high 
deflection  angles  and  was  proven  inappropriate.  In  high  speed 
theory  based  on  Eckert's  reference  temperature,  low  tempera¬ 
ture  viscosity  shov/s  no  divergence.  High  speed  theory  is 
valid  according  to  test  conditions  and  correlation  vrith  numer¬ 
ical  results.  The  numerical  solution  and  Eckert  theory  yield 
the  isothermal  wall  solution  used  for  comparison, 

numerical  Results.  Numerical  calculation  for  the  noniso- 
thermal  wall  shov/ed  a  divergence  from  the  isothermal  case  ac¬ 
counting  for  28  percent  of  the  discrepancy  cited  at  a  deflec¬ 
tion  angle  of  oC  =  3  degrees  ,  Originally,  this  calculation 
also  accounted  for  100  percent  of  the  discrepancy  at  oC=  14 
degrees  .  Hov/ever,  according  to  the  study  accomplished  by 
Woo  (Ref  9),  the  experimental  data  consistently  yielded  a 
ratio  of  heat  transfer  coefficient  to  reference  of  = 

2,4  •  With  this  discrepancy  value,  the  nonisothermal  wall 

accounted  for  41  percent  of  the  deviation.  The  step  size 
in  the  above  cases  was  derived  directly  from  the  thin  skin 
and  HR3I  test  data.  The  step  size  was  varied  to  model  pos¬ 
sible  wall  temperature  conditions.  The  final  wall  temperature 
was  varied  by  doubling  and  quadrupling  the  initial  step  size 
and  by  doubling  the  initial  wall  temperature  for  all  deflec¬ 
tion  angles  considered.  The  result  of  increasing  magnitude 


of  the  step  is  a  significant  decrease  in  heat  transfer  coef¬ 
ficient,  and  increase  in  the  downstream  distance  affected  by 
the  wall  temperature  discontinuity. 

The  initial  wall  temperature  was  varied  to  model  the  lower 
surface  temperature  expected  on  the  stainless  steel  wedge. 

The  initial  v/all  temperature  assumed  the  value  of  the  initial 
thin  skin  temperature,  with  the  final  wall  temperature  taken 
from  the  appropriate  HRSI  data.  These  runs  v/ere  expected  to 
best  correlate  the  actual  wall  temperature  discontinuity  in 
the  wind  tunnel  tests.  The  result  shov/ed  the  nonisothermal 
wall  effect  accounted  for  40  percent  of  the  discrepancy  at 
0(  =  3  degrees  ,  and  60  percent  at  oC=  14  degrees  .  Ac¬ 
cordingly,  the  wall  temi)erature  step  v/ould  have  to  be  twice 
its  expected  value  to  account  for  the  entire  discrepancy. 

The  location  of  the  step  was  varied  to  investigate  the 
effect  of  flov.’  development  on  the  heat  transfer  coefficient. 
The  divergence  of  the  heat  transfer  coefficient  from  the  iso¬ 
thermal  case  grew  with  increasing  distance  dovmstream,  Re- 
s’‘onse  also  seemed  slov/er  for  the  aft  step  locations. 

The  results  of  the  numerical  calculations  are  valid  as 
long  as  the  boundary  layer  assumptions  are  valid,  Ho'.^ever,  in 
the  area  of  the  v;all  temperature  discontinuity,  the  boundary 
layer  solution  is  not  valid  due  to  axial  diffusion  and  three 
dimensional  effects.  In  the  area  of  the  discontinuity,  the 
numerical  results  are  questioned.  Away  from  this  disconti¬ 
nuity,  the  boundary  layer  solution  should  be  reliable.  The 
effect  of  diffusion  near  the  discontinuity  nay  be  a  delay  in 
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recovery  of  the  heat  transfer  coefficient  to  a  point  further 
dovmstream.  In  addition,  three  dimensional  effects  at  the 
surface  around  the  discontinuity  may  keep  the  tile  surface 
temperature  closer  to  the  original  wall  temperature.  The  ef¬ 
fective  step  may  be  delayed  such  that  the  location  of  the 
thermocouple  may  be  in  the  region  of  steep  heat  transfer  co¬ 
efficient  recovery.  The  combined  effects  of  these  character¬ 
istics  would  be  a  significantly  lov/er  heat  transfer  coefficient 
revading  at  the  therraocouple  location,  possibly  still  account¬ 
ing  for  the  entire  discrepancy. 

Analytical  Results.  Analytical  solution  of  the  boundary 
layer  with  the  nonisothermal  v/all  condition  was  accomplished 
by  wall  temperature  superposition.  The  superposition  method 
revealed  similar  effects  of  the  nonisothermal  wall  to  those 
shoTO  by  the  numerical  solution.  However,  some  magnitude 
discrepancy  was  shown  between  the  numerical  and  analytical 
methods.  Again,  the  analytical  methods  assume  boundary  layer 
conditions  and  nay  not  be  valid  near  the  v/all  temperature  dis¬ 
continuity. 

Discussion  of  Results.  The  results  of  this  study  with 
respect  to  the  original  problem  posed  are  summarized.  The 
nonisothermal  v;all  effect  does  account  for  much  of  the 
perinental  data  discrepancy,  with  as  much  as  half  of  the  dis¬ 
crepancy  accounted  for.  Since  exact  surface  temperature  dis¬ 
tribution  is  unknovm  on  the  KP3I  tile,  the  wall  temperature 
step  is  still  not  completely  'lefined.  The  lov;est  possible 
temperature  on  the  stainless  steel  wedge  is  the  thin  skin 
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initial  temperature.  The  wedge  surface  is  expected  to  have 
been  very  close  to  this  temperature  throughout  the  run.  The 
step  size  was  probably  very  close  to  that  used  in  the  initial 
v;all  temperature  variation  calculations.  This  result  showed 
the  nonisothermal  wall  accounting  for  roughly  half  of  the 
original  discrepancy.  The  ITR3I  temperature  may  have  been 
somewhat  higher  on  the  surface  than  the  data  indicated.  The 
possible  delay  in  heat  transfer  coefficient  response  due  to 
axial  diffusion  and  due  to  three  dimensional  effects  at  the 
surface  may  be  an  additional  effect  not  accounted  for.  Since 
the  thermocouple  location  is  near  the  v;all  temperature  step 
location,  the  tem'^erature  reading  may  have  been  in  the  region 
of  the  heat  transfer  coefficient  recovery,  A  much  lov/er  value 
of  this  coefficient  could  result.  -lo'wever,  the  final  conclu¬ 
sion  of  these  findings  is  that  only  one-half  of  the  discrep¬ 
ancy  may  actually  be  attributed  to  the  nonisothermal  effect, 
pending  further  investigation.  The  remaining  port'on  may  be 
due  to  the  diffusion  effect,  three  dimensional  effects,  in¬ 
accuracy  in  thermocouple  data,  fluctuations  in  boundary  layer 
and  surface  heating  ahead  of  the  tile,  and/or  inaccuracy  in 
data  reduction,  'hile  the  radiation  heat  transfer  v?as  also 
cited  as  a  possible  cause  earlier,  the  IIPSI  temperature  should 
not  vary  greatly  from  the  thermocouple  data.  As  a  result, 
radiation  analysis  fnould  serve  to  verify  the  test  data,  at 
be  st , 

Analicatl ~ns.  The  results  of  this  study  m^y  be  applied 
beyond  the  scope  of  the  original  problem.  The  effect  of  the 
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nonisothermal  wall  is  present  in  a  number  of  applications, 
particularly  in  high  speed  flow,  on  any  high  speed  vehicle, 
discontinuity  in  surface  material  results  in  temperature  dis¬ 
continuity,  This  is  especially  true  of  the  Space  shuttle 
Orbiter,  Several  locations  are  marked  by  a  change  in  surface 
material  on  this  vehicle.  At  the  nose cap,  the  surface  is 
insulated  vrith  Reusable  Carbon-Carbon  (RCC)  material,  follov/ed 
by  HRjI  surface  insulation.  The  Carbon-Carbon  acts  as  a  heat 
sink  v.'hile  the  ERSI  is  an  insulator.  As  such,  the  initial 
temperature  should  be  lower  than  the  HP.3I  surface  temperature, 
and  a  wall  temperature  step  is  present.  The  data  presented 
here  qualitively  describes  the  fluctuation  of  heating  due  to 
this  discontinuity,  depending  on  the  relative  magnitude  of  the 
RCC  and  Tf.'jl  surface  temperatures,  similar  surface  material 
changes  occur  at  other  locations.  The  lov/er  surface  points 
of  attachment  to  the  External  Tank  ("JT)  are  marked  by  an  IB3I, 
R.CC,  stainless  steel,  RCC,  HR3I  sequence  of  materials.  The 
heating  characteristics  near  these  points  should  reflect  a 
compounded  nonisothermal  \ra.ll  effect.  On  the  0"3  pod,  a  low- 
temperature  Reusable  Surface  Insulator  (LRSI)  to  ?R3I  material 
change  occurs.  Again,  a  wall  temperature  step  occurs,  v:ith 
the  associated  effects.  Of  course,  the  flov?  is  further  con- 
rlicated  by  changing  boundar5’^  layer  edge  conditions  over  the 
vehicle,  but  the  demonstrated  nonisothermal  v'all  effect  v;ill 
greatly  influence  the  surface  heating  at  many  locations  on 
+he  vehicle. 
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Recommendations 

As  discussed  above,  numerous  questions  remain  unanswered. 
The  temperature  distribution  is  unlcnovm  on  the  surface  of  the 
stainless  steel  v/edge.  In  order  to  completely  define  the 
v/all  temperature  and  investigate  the  full  effects,  the  sur¬ 
face  temperature  ahead  of  the  test  article  location  should  be 
investigated.  This  may  be  accomplished  through  experimental 
tests  at  all  deflections  using  thermocouple  measurements,  or 
through  a  thorough  analytical  model  of  the  material  conduc¬ 
tive  and  radiative  properties  based  on  flov;  conditions.  The 
result  of  this  investigation  v/ould  be  better  knov:ledge  of 
the  temperature  distribution  and  subsequent  effects  of  the 
wall  temperature  variation  on  heating.  The  HH3I  surface  tem¬ 
perature  should  also  be  verified,  since  the  test  thermocouple 
data  was  taken  beneath  the  radiative  coating.  In  addition, 
knov/ledge  of  the  surface  temperature  of  the  HRdl  tile  near 
the  surface  material  change  would  indicate  the  location  of 
the  effective  step  and  help  to  further  model  the  heat  transfer 
coefficient  variation.  AC'in,  thermocouple  test  data  could 
be  used. 

The  wall  temperature  of  the  v/ind  tunnel  during  the  experi¬ 
mental  runs  v;ould  have  served  to  validate  the  TIR3I  surface 
temperature  distribution.  V/ithout  knovrledge  of  the  tunnel 
wall  temperature,  equilibrium  heating  could  not  be  fully  de¬ 
termined,  This  prevented  knov/ledge  of  surface  tempera¬ 

ture,  especially  at  the  higher  deflection  angles.  This  should 
be  investigated  by  experimental  observation. 
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The  boundary  layer  assumptions  were  considered  invalid  at 
the  discontinuity.  As  a  result,  the  results  near  the  discon¬ 
tinuity  may  not  be  reliable,  A  complete  Kavier-Stokes  solu¬ 
tion  of  the  boundary  layer  should  be  accomplished  to  determine 
the  limits  of  validity  of  the  boundary  layer  solution,  to 
verify  the  nunorical  results  at  all  locations,  and  to  investi¬ 
gate  further  the  com’jlete  nonisothermal  v/all  effect. 

Possible  errors  cited  include  thermocouple  inaccuracy, 
data  reduction  in-accuracy,  and  conduction  in  the  tile.  The 
reliability  of  thermocouple  data  in  high  speed  testing  is 
inherent  to  the  accuracy  and  validity  of  test  results  in  all 
applications,  prrors  caused  by  mounting,  location,  and  con¬ 
duction  shovild  be  investigated.  Conduction  in  the  tile  ma¬ 
terial  is  perhaps  the  most  important  property  of  the  insulator 
since  it  determines  the  ca:ability  to  protect  the  Orblter  dur¬ 
ing  reentry.  All  conductive  effects  must  be  understood,  par¬ 
ticularly  those  due  to  irregularities  in  the  coating,  in  mount¬ 
ing,  r.nd  in  construction  of  the  tile  itself.  Purther  investi¬ 
gation  is  needed  in  this  area.  The  study  of  'oo  (Pef  9)  has 
served  to  verify  the  accuracy  of  the  data  reduction  technique. 
Validation  of  the  technique  for  knovm  results  may  further 
verify  its  reliability. 

Lastly,  many  apli  cat  ions  of  the  nonisothermal  v;all  exist, 

A  parametric  study  of  v;all  temperature  variation  effects  on 
heating  should  be  conducted  v;ith  the  emphasis  on  applications. 
This  '..'ould  encompass  a  broad  range  of  flov:  conditions,  step 
locations,  step  sizes,  and  temperature  distributions.  The 


result  of  this  type  of  study  would  be  enhancement  of  applica¬ 
bility  of  results  to  varied  conditions  and  to  a  broader  range 
of  fields.  The  nonisothermal  wall  effect  on  heating  could 
thus  be  better  understood,  recognized,  and  predicted  in  future 
applications. 


Table  I 

HR3I  Tile  Run  Description 

Angle  at 

Pile  Injection  Description 


1 

14  degrees  . 

.  .hold 

.  .  9  deg 

• 

•  •  holcl  •  • 

14  deg 

2 

14  degrees  . 

.  .hold 

,  .  9  deg 

• 

*  •  hoXd.  •  • 

14  deg 

3 

14  degrees  . 

.  .hold 

.  .  9  deg. 

• 

.no  hold.  , 

14  deg 

4 

14  degrees  , 

.  .hold 

.  .  3  deg. 

• 

•  • hold  •  • 

14  deg 

5 

14  degrees  . 

.  . hold 

.  .  3  deg. 

• 

.no  hold.  , 

14  deg 

6 

10 

degrees  , 

.  3  deg 

.  .  14  deg 

* 

.  3  deg  .  . 

14  deg 

7 

14  degrees  . 

.  .hold 

.  .  3  deg 

• 

•  •hold.  •  • 

14  deg 

8 

14 

degrees  . 

.  .hold 

.  .  9  deg. 

• 

.no  hold,  . 

14  deg 

9 

Fitch  error, 

file  scratched. 

10 

14 

degrees  . 

.  .hold 

.  .  9  deg. 

• 

•  • hold  •  • 

14  deg 

11 

10 

degrees  , 

.  3  deg 

.  ,  hold  . 

• 

•  • hold  •  • 

10  deg 

12 

14 

degrees  , 

.  .hold 

.  ,  hold  . 

• 

.  .hold  .  . 

hold 

13 

14 

degrees  . 

.  .hold 

,  .  3  deg. 

• 

•  • hold  •  • 

14  deg 

14 

10 

degrees  , 

.  3  deg 

.  ,14  deg. 

• 

.  3  deg  .  . 

14  deg 

15 

10 

degrees  , 

.  .hold 

.  .  shock  generator  .  . 

hold 

16 

10 

degrees  , 

.  3  deg 

.  .  hold  . 

shock  generator  .  . 

17 

14 

degrees  . 

.  .hold 

.  ,  shock  generator  ,  .  hold 

18 

14 

degrees  . 

.  .hold 

.  .  3  deg. 

• 

.  .hold  ,  , 

14  deg 

19 

14 

degrees  . 

•  .hold 

.  .  hold  . 

• 

•  • hold  •  • 

•  •  • 

20 

10 

degrees  . 

.  . hold 

,  ,  hold  , 

• 

•  • hold  •  • 

•  •  • 

21 

10 

degrees  . 

.  3  deg 

.  ,  hold  . 

• 

•  • hold  •  • 

•  •  • 

22 

14 

degrees  . 

.  .hold 

,  ,  9  deg. 

• 

.  .hold  ,  , 

14  deg 

23 

14 

degrees  . 

.  .hold 

,  ,  9  deg. 

• 

.  ,14  deg  . 

•  •  • 
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Table  II 

FR3I  Tile  Run  Description 


Pile 

Angle  at 
Injection 

Description 

27 

14 

degrees 

• 

• 

.  hold  . 

• 

• 

hold  ...  hold  . 

•  • 

28 

10 

degrees 

• 

• 

.  hold  . 

• 

• 

hold  .  ,  .  hold  . 

•  • 

29 

10 

degrees 

• 

• 

.  3  deg 

• 

• 

hold  ...  hold  . 

•  • 

30 

10 

degrees 

• 

• 

.  hold  , 

• 

• 

hold  .  .  .  hold  . 

•  • 

31 

14 

degrees 

• 

• 

,  hold  . 

• 

9  deg  .  .  hold  . 

14  deg 

32 

14 

degrees 

• 

• 

,  hold  . 

• 

♦ 

9  deg  .  .  14  deg 

•  • 

33 

14 

degrees 

• 

• 

.  hold  , 

• 

• 

5  deg  .  .  hold  . 

14  deg 

34 

14 

degrees 

• 

• 

no  hold 

• 

• 

3  deg  .  .  14  deg 

hold 

35 

14 

degrees 

• 

• 

no  hold 

• 

• 

3  deg  .  .  hold  . 

11.6  deg 

36 

10 

degrees 

• 

• 

.  3  deg 

# 

« 

hold  .  shock  generator  .  • 

37 

10 

degrees 

• 

• 

.  hold  . 

• 

• 

shock  generator  , 

#  » 

38 

14 

degrees 

• 

• 

,  hold  . 

• 

• 

shock  generator  . 

•  • 

39 

14 

degrees 

« 

• 

.  hold  . 

# 

• 

3  deg  .  .  hold  . 

14  deg 

40 

14 

degrees 

• 

• 

,  hold  . 

• 

• 

3  deg  .  no  hold 

14  deg 

41 

14 

degrees 

• 

• 

,  hold  , 

• 

• 

hold  .  .  .  hold  . 

•  • 

42 

10 

degrees 

• 

• 

,  hold  , 

• 

• 

shock  generator  . 

•  • 

43 

10 

degrees 

• 

• 

,  hold  , 

• 

• 

shock  generator  . 

•  • 

44 

10 

degrees 

• 

• 

,  hold  . 

• 

• 

shock  generator  . 

•  • 

45 

14 

degrees 

• 

.  hold  . 

• 

• 

shock  generator  . 

•  • 

44 


Table  m 

Thin  Skin  Run  Description 


Angle  at 

Pile  Injection  Description 


46 

14 

degrees 

• 

• 

• 

hold  . 

• 

.  hold  .  ,  ,  hold 

• 

• 

• 

47 

14 

degrees 

• 

• 

• 

hold  . 

• 

shock  generator 

• 

• 

• 

• 

48 

14 

degrees 

• 

• 

• 

hold  . 

• 

,  hold  ...  hold 

m 

« 

• 

49 

14 

degrees 

• 

• 

• 

hold  . 

• 

shock  generator 

• 

• 

• 

50 

14 

degrees 

• 

• 

• 

hold  , 

• 

shock  generator 

• 

• 

• 

• 

51 

14 

degrees 

• 

• 

• 

hold  . 

• 

shock  generator 

• 

• 

• 

• 

52 

14 

degrees 

• 

• 

• 

hold  , 

• 

shock  generator 

• 

• 

• 

• 

53 

14 

degrees 

• 

« 

• 

hold  . 

« 

shock  generator 

• 

• 

• 

• 

54 

10 

degrees 

• 

• 

• 

hold  , 

• 

shock  generator 

• 

• 

• 

• 

55 

10 

degrees 

• 

• 

• 

hold  . 

• 

shock  generator 

• 

• 

• 

• 

56 

10 

degrees 

« 

• 

• 

hold  . 

• 

shock  generator 

• 

• 

» 

• 

57 

10 

degrees 

• 

• 

0 

hold  . 

• 

.  hold  .  .  .  hold 

• 

• 

• 

58 

10 

degrees 

• 

• 

• 

hold  . 

• 

.  3  deg  .  .  hold 

• 

• 

• 

59 

10 

degrees 

• 

• 

• 

3  deg 

• 

,  hold  .  shock  generator 
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7” 


.16"  SIP  .125"  Aluminum 


Figure  1.  HR3I  Test  Apparatus 
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OL  (deg) 


Figure  5.  Heating  Estimate  for  Wind  Tunnel  Thermocouple 
Data  (Transient  Maneuver)  (Ref  1) 


Figure  ?•  Thermocouple  Data  for  rosi  Test  Runs 


Figure  8,  Thermocouple  Data  for  Thin  skin  Test  Runs 
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Figure  9,  'rfind  Tunnel  Thernocouple  Data  Distribution 
{Cx=  3  degrees) 


Figure  10.  w^d  Tunnel  Thermocouple  Data  Distribution 

=  10  degrees) 


Figure  1 1 


’,7ind  Tunnel  Thermocouple  Data  Distribution 
{OC  =  H  degrees) 


1 


O  Run  58,0{=  3  degrees 
^  Run  57,Of=  10  degrees 
□  Run  A3fOC=  14  degrees 


0  L 
7 


Figure  12 
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X  (in.) 

,  Experimental  Heat  Transfer  Coefficient  Distri 
bution  from  Thin  Skin  Wind  Tunnel  Tests 


T 


(deg  R) 


fl  (xIO"*^  Ib^sec/ft^) 


Figure  15.  Comparison  of  Sutherland's  Law  and  Low 
Temperature  Calculations  of  Viscosity 


OC  (deg) 

Figure  16,  comparison  of  isothermal  ’./all  Heat  Transfer 
coefficient  calculations  (x/L  =  0.575) 


O  Numerical  Solution 
O  High  Speed  Eckert  Theory,  Eq  (15) 
Eckert  Theory,  Eq  (10) 


h 

(x10“'^ 

Btu 

ft*s.R 


x/l 


Figure  21, 


Streamwise 

Coefficient 


Distribution  of  Heat  Transfei 
(Of =  14  degrees) 
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Figure  28,  Comparison  of  Heat  Transfer  Coefficient  with 
Varying  Step  Size  (x/l  =  0,575) 


Figure  29.  Comparison  of  Streamv;ise  Distribution  of 
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Figure  55.  Comparison  of  Nonisothernal  Wall  Calcula¬ 
tions  of  Heat  Transfer  Coefficient 
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Biblio/;raphy 


1.  Hodge,  J,  K. ,  P.  V/.  Phillips,  and  D.  R.  Audley.  "Plight 
Testing  a  Manned  lifting  Reentry  Vehicle  (Space  Shuttle) 
for  Aerothermodynamic  Performance,”  AIAA  Paper  AIAA-81- 
2421,  AIAA  1st  plight  Testing  Conference,  November  1981. 

2.  Eckert,  E.  R.  G.  "Survey  of  Boundary  Layer  Heat  Trans¬ 
fer  at  High  Velocities  and  High  Temperatures,”  V/ADC 
Technical  Report  59-624,  WADC,  April  I96O. 

3.  Hayes,  j.  R.  and  R.  D.  Neumann,  "Trends  in  the  Space 
Shuttle  Aerothermodynamic  Data  Base,"  AIAA  Paper  AIAA- 
81 -0292,  AIAA  19th  Aerospace  sciences  Meeting,  January 

1981. 

4.  Kays,  N.  M.  Convective  Heat  and  Mass  Transfer.  New 
York ;  McGrav;-Hill  Book  Company,  1 966, 

5.  Piore,  A.  V/.  "Viscosity  of  Air,"  journal  of  space¬ 
craft  and  Rockets.  3  (5):756-758  (May  1966). 

6.  Keyes,  p.  G.  "A  Summary  of  Viscosity  and  Heat  Conduc¬ 
tion  Data  for  He,  A,  ^2$  N2»  CO,  C02»  H2O,  and  Air," 

Transactions  of  ASMS.  73:589-596  (1951). 

7«  Hayes,  J.  R,  Data  transmittal  of  magnetic  tape  and  thin 
skin  data  from  shuttle  test.  Air  porce  Nright  Aeronautical 
Laboratories,  Nright-patterson  APB,  13  April  1981, 

8,  Holman,  j,  p.  Heat  Transfer  (pifth  edition).  New  York: 
McGraw-Hill  Book  Company,  1981, 

9.  Noo,  Y.  K.  Transient  Heat-Transfer  Measurement  Tech¬ 
nique  in  Vina  Tunnel  and  Data  Analysis  Technique  Using 
System  Identification  Theory.  MS  Thesis.  Nright- 
Patterson  A?B,  Ohio:  Air  porce  Institute  of  Technology, 
December  1982, 

10.  Hankey,  'Vilbur,  Lecture  notes  distributed  in  AE  6,50, 
Reentry  Aerodynamics,  School  of  Engineering,  Air  porce 
Institute  of  Technology,  v/right-patterson  APB,  April 

1982. 

11,  Cohen,  C.  B.  and  E.  Reshotko.  "similar  Solutions  for 
the  Compressible  Laminar  Boundary  Layer  with  Heat  Trans¬ 
fer  and  pressure  Gradient,"  NACA  Report  1293,  1956, 


79 


12,  Nachtsheim,  P,  R,  and  P,  Swigert,  "Satisfaction  of 
Asymptotic  Boundary  Conditions  in  Numerical  Solution  of 
Systems  of  Nonlinear  Equations  of  Boundary  Layer  Type," 
NASA  Report  TN  D-5004,  October  1965# 

13.  Beauregard,  A.  J.  An  Analytical  Study  of  the  Effects 
of  Nass  Transfer  on  a  Gomoressible  Turbulent  Boundary 
Layer.  IIS  Thesis.  V/righi-patterson  ARB,  Ohio:  Air 
Force  Institute  of  Technology,  December  1976* 


Appendix  A 


Similarity  Derivation  and  Applicability 

Since  much  of  the  focus  of  this  project  v;as  based  on  a 
model  of  the  boundary  layer  flov;,  the  application  of  simi¬ 
larity  equations  was  a  method  worth  consideration.  No  ref¬ 
erence  was  found  yielding  the  derivation  of  similarity  equa¬ 
tions  for  high  speed  flow  v/ith  Prandtl  Nunber  unequal  to  unity. 
For  future  consideration,  this  derivation  is  included  here, 
and  v;ill  parallel  the  development  used  by  Hahkey  (Ref  10), 
Starting  \/ith  the  tv;o-dimensional  boundary  layer  equa¬ 
tions, 


=  0  (A-1  ) 

/>uu^  +  fvuy  =  (^Uy)y  .  (A-2) 

*  ^THy  =  aKV  [( '  -W ) Wy]  y  (  A-5  ) 

and  defining  the  stream  function  1^  such  that  p\x  =’]^y  ^•nd 
pv  =  ,  the  boundary  layer  equations  may  be  written  as 

>xy->yx  =  ° 

T^y^x  “  *^x^y  ~  ^e^e^e  (A-5) 

+  'S'V  C  -  *^>''Ty]y  (A-6) 


where  gq  (A-4)  is  identically  satisfied  as  long  as  the  stream 
function  is  continuous  v;ith  resncct  to  the  independent  vari- 
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ables  X  and  y.  If  the  stream  function  is  assumed  separable 
v/ith  respect  to  the  transformed  variables  ^  and>J  ,  it  may  be 
written  as 


y  =  Sit)  f('n) 


(A-7) 


where  g  =  v/?T  t  and  the  transformed  variables  and  are 
related  to  the  independent  variables  x  and  y  as 


df  =  j(x)dx 


CA-8) 


dTJ  =  h(,x)pdy 


(A-9) 


with  3  =  and  h  =  ^g/s  •  By  chain  rule,  differ¬ 

entiation  v;ith  respect  to  the  variables  x  and  y  may  now  be 
defined  as 


3  If  In 


5y  = 


(A-10) 

(A-11) 


Making  use  of  the  above  transformations  and  definitions,  Eq 
(A-5)  is  reduced  to 

.2 _  ^  .  .2„2. 


Defining  0  =  ^  allows  further  simplification  to 
>^/^e 


hpUeCTv^Uj  =  h2f(cu^)^+  |>gUgUg^ 

Dividing  through  by  h^  and  reducing  yields 


(A-13) 


^  (A-u) 
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i’jsswtKSf  c 


I 


letting  f*  =  ^-^1=  u/Ug  and  S  =  (h/h^)  -  1  gives 

(Cf*)'  =  g[gf  ^  ^  (A-15) 


Further, 


Hi  =f.  .H_!ki  +  lii.H,u 

“e  “e  *  “e  “e  “e  ?  H 


(A-i6) 


giving  the  final  form 


(cf'')‘  +  gf'V'f  =  ^  (A-1T) 


similarly,  the  development  of  3<1  (A-6)  is  found  to  yield 
"^y^x  "  ^^y  “  ■*■  Pr^  ^yJ  y  ' 


(A-18) 


which,  applying  the  transformation  to  variables  %  nfy  gives 


;[ef -Hj  -  IjhJ  =f^  (uu^+  ^  TV  ^ 


(A-19) 


With  the  definition  f^T  s,  the  total  enthalpy  H  may  be  used 


to  transform  the  temperature  T; 


H  =  (S  +  1)Hg 
H  -  u^/2 


uu^) 


(A-20) 


(A-21 ) 


(A-22) 


Differentiating  the  total  enthalpy  and  velocity  as 


H,y^=  HeS* 


=  Hg^(S  +  1)  +  HgS^ 


(A-23) 


(A-24) 


^11=  V” 


(A-25) 


These  values  can  be  substituted  into  Eq  (A-19)  to  give 

g^f’[Hej(S  +  1)  +  HgSj]  -  gHgS'tj  = 


(A-26) 


PM^ 

Again  using  the  definition  of  C  =  and  realizing  that 

total  enthalpy  is  constant,  Eq  (A’~2S)  becomes 

2 

(CS'^+  PrgS't^  =  Pr[g^f‘S^  -  ^  (^^)(Cf 'f" )J(A-27) 

Applying  the  similarity  conditions  that  the  value  C  is  con¬ 
stant  and  equal  to  unity  and  that  the  functions  f  and  S  are 
independent  of  strearawise  location,  (i,e.  Sj  =  0  ,  1^  = 

f/g  )»  Sqs  (A-17)  and  (A-27)  are  further  reduced  to 


f  •  •  +  ff«  =/3(f»^  -  S  -  1 ) 

S"  +  PrfS*  =  (1  -  Pr)f-^^7y%l(f'f"  ’  +  f’^) 


(A-28) 


(A-29) 


where  j6  =  g^  ^  and  f ^  ^  (f -  S  -  1 ) 

Note  that  these  expressions  agree  with  those  given  by  Cohen 
and  Reshotko  (Ref  11).  These  equations  are  valid  so  long  as 
the  similarity  conditions  are  met. 

For  the  conditions  of  the  wind  tunnel  tests  run,  the  in¬ 
sulating  articles  used  were  assumed  to  be  at  thermal  equili¬ 
brium  between  boundary  layer  heating  and  radiation.  For  radi¬ 
ation 


q  =tCrT5 


(A-30) 


assuming  zero  sink  temperature.  The  heating  of  the  wall  is 
given  as 


At  the  wall,  velocity  is  zero  and  H  =  c^T  or  T  =  e/o^ 
substituting  and  equating  the  above  expressions  yields 


(A-31) 


p  Jy=o 


(A-32) 


Defining  the  function  S  as  before  and  using  the  differentia¬ 
tion  of  Sqs  (A-10)  and  (A-11),  the  equilibrium  condition  is 
expressed  as 


'“e/’e'fe'=p 


(S  +  1) 


(A-33) 


in  this  form,  g  is  a  function  of  the  streamwise  location  de¬ 
noted  by  i.  As  a  result,  the  function  3  is  no  longer  inde¬ 
pendent  of  "f,  as  necessary  by  the  similarity  assumptions. 

In  this  case  similarity  is  no  longer  valid,  and  must  be  dis¬ 
carded  as  a  solution  technique. 

For  those  cases  v/here  similarity  is  valid,  solution  of  the 
given  equations  may  be  accomplished  according  to  the  tech¬ 
nique  given  by  Kachtsheim  anci  fSwigert  ^Ref  12), 
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Appendix  B 


Wumerical  Commutation 

All  numerical  calculations  were  accomplished  using  the 
boundary  layer  program  alluded  to  earlier.  Access  to  this 
program  was  gained  from  nr,  George  Havener,  The  program  used 
was  identical  to  that  found  in  Ref  13  given  the  folowing 
modifications : 

1)  omit  lines  200-220, 

2)  omit  line  260, 

3;  substitute  for  line  420  with 

READ  (5,8002)  G,  PR,  XMINP,  R3Y,  TA 

4)  omit  lines  480-610, 

5)  omit  lines  5120-5130, 

6)  omit  lines  5260-5310, 

7)  substitute  between  lines  5250  and  5320  with 

D3(1)  =  0. 

8)  substitute  between  lines  7480  and  7490  v;ith 

TTR  =  (TA  +  112.)/(TA*TREP  +  112,) 

9)  omit  lines  7990  and  8020, 

10)  substitute  betv/een  lines  7980  and  8030  with 

EP(K)  =  EP(IT)*(1.-EXP(-.412*((S-BTRX)/X1AII)**2,)) 

11)  omit  lines  11950-12450. 

The  program  itself,  called  ITRACT,  computes  boundary  layer 
flov^  characteristics  for  laminar  and/or  turbulent  boundary 
layers,  over  flat  plates  or  axisymmetric  bodies.  The  com¬ 
putation  was  accomplished  through  trsmsforraation  of  the  non¬ 
linear  parabolic  boundary  layer  eciuations  into  a  rectangular 
grid,  and  subsequent  solution  of  this  system  of  equations 
by  linear  finite  differencing,  using  a  three-point  scheme 
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with  a  tridiagonal  matrix  inversion  routine.  The  solution 
of  the  boundary  layer  was  thus  calculated  at  incremental 
streamwise  locations.  A  complete  discussion  of  the  technique 
used  in  this  boundary  layer  program  is  given  by  Beauregard 
(Ref  15). 

The  characteristics  input  to  the  program  itself  included 
the  specific  heat  ratio  JT,  Prandtl  Number  Pr,  boundary  layer 
edge  Mach  Number  M. ,  edge  temperature  T.,  edge  Reynolds  Number 
based  on  total  plate  length  Rg2»  wall  temperature 

to  stagnation  temperature  ratio,  flov;  transition  location, 
viscosity  lav;  selection  parameter,  and  calculation  parameters. 

The  specific  heat  ratio  and  prandtl  Number  were  assumed  con¬ 
stant  for  all  cases  at  y=:  1,4  and  Pr  =  0.72  .  Boundary 

layer  edge  properties  of  Mach  Number,  static  temperature,  and 
Reynolds  Number  v;ere  input  appropriately  for  the  case  considered. 
V/all  temperature  to  stagnation  temperature  ratio  was  derived 
from  the  v;all  temperature  given  by  thin  skin  data,  assuming 
a  stagnation  temperature  of  T^  =  2000  R  for  all  cases. 
Sutherland's  Lav;  for  viscosity  was  used  in  all  cases,  except 
as  discussed  in  consideration  of  the  alternate  low  temperature 
viscosity  relationship.  The  flov;  transition  location  was 
arbitrarily  chosen  as  unity  (end  of  the  plate)  since  flov;  was 
assumed  to  be  laminar  in  all  cases.  This  assumption  was  made 
for  several  reasons.  No  indication  of  transition  over  the 
test  apparatus  was  found,  and  no  location  for  transition 
could  be  given.  The  freestream  Reynolds  Number  was  found 
to  be  approximately  one  million,  with  the  location  of  interest 
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at  approximately  one  half  the  plate  length  giving  a  Re3rnolds 
ITumber  locally  of  roughly  one-half  million.  In  addition, 
much  care  was  taken  in  installing  the  test  articles  to  pre¬ 
vent  a  boundary  layer  trip.  Since  heating  data  for  the  test 
runs  made  give  relatively  smooth,  continuous  heating,  no 
transition  is  indicated.  Lastly,  experience  with  the  wind 
tunnel  used  indicates  that  turbulent  flow  is  unlikely  and 
difficult  to  attain  at  operational  conditions. 

Output  data  included  boundary  layer  velocity  profile, 
thickness,  momentum  thickness,  displacement  thickness,  local 
Reynolds  ijumber,  and  local  Stanton  Number,  prom  the  local 
Stanton  Number,  the  local  heat  transfer  coefficient  was  de¬ 
rived. 

In  order  to  estimate  the  effect  of  the  wall  temperature 
step  on  the  heat  transfer  coefficient,  the  program  had  to  be 
modified  to  accommodate  a  step  change  in  the  v/^all  to  stagna¬ 
tion  temperature  ratio.  This  modification  was  accomplished 
by  inserting  the  following  lines  between  lines  1480  and  1500 
in  the  program  listed  in  Ref  12: 

IP  (S.  LT  .  0.467)  SO  TO  500 
BO  =  .5500 
T’:  =  BO  *  T10 
500  DX2D3  =  DX1D3 

The  value  of  the  v;all  to  stagnation  temperature  ratio,  denoted 
BO  above,  is  based  on  the  total  wall  temperature  following 
the  step,  given  by  the  HRSI  test  data  cited.  This  value 
varies  depending  on  the  stagnation  temperature  and  final  wall 
temperature  desired. 
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Sstimation  of  boundary  layer  characteristics  based  on  the 
low  temperature  viscosity  calculation  required  some  additional 
modification  of  the  program.  The  following  changes  account 
for  the  alternate  viscosity  relationship: 
omit  line  880, 

2)  substitute  the  follov;ing  lines  betv/een  lines  870  and 
890: 

SU  =  220./10.**(9./TA) 

TR  =  SU/TA 

TP.EP  =  (G-1.)*XMINI^*2. 

TC  =  TR/TREP 

5)  substitute  the  following  for  lines  900  to  910: 

EPS  =  (((1.+TR)*TREP**1.5)/((TREP+TR)*RSY))**0.5 

4)  substitute  for  line  930  with 

102  COrTIMUE 

5)  omit  line  1650, 

6}  substitute  for  line  1670  with 

676  xiTus  =  VI  Sir? 

7)  omit  line  1300. 

Note  that  the  modifications  for  the  lov;  temperature  viscosity 
calculation  assume  that  this  method  is  replacing  the  Suther¬ 
land*  s  Law  calculation  and  replaces  statements  accordingly. 

The  same  viscosity  calculation  parameter  is  used  in  the  in¬ 
put  data.  Laminar  flow  is  also  assumed,  and  additional  changes 
may  be  required  for  turbulent  flow  applications. 

Lastly,  a  check  was  made  for  numerical  solution  inac¬ 
curacy  due  to  increment  size.  This  check  was  made  for  the 
nonisothermal  case  of  OC  =  3  degrees  .  The  increment  size, 
originally  set  at  Ax  =  0.0001  ,  was  doubled  to  Ax  =  0.0002  , 

and  halved  to  Ax  =  0.00005  to  verify  accuracy.  The  result 
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for  this  check  showed  an  error  of  five  percent  above  the 
value  used  both  before  the  step  and  at  the  HRSI  thermocouple 
location  for  the  increment  Ax.  =  0.0002  •  When  the  increment 

was  halved,  or  Ax.  =  0,00005  *  an  error  of  two  percent  be¬ 
low  the  calculated  value  was  found  between  this  and  the  ori¬ 
ginal  case  of  Ax  =  0.0001  ,  for  both  locations,  before  and 

after  the  v/all  temperature  step.  It  would  appear  that  accu¬ 
racy  was  still  increasing  with  decreasing  increment  size  at 
the  values  used. 
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(from  Item  20)^using  a  numerical  boundary  layer  program,  supplement¬ 
ed  by  wall  temperature  superposition  analysis.  The  program  was 
shown  to  cive  results  consistent  with  Bckert  high  speed,  flat  plate 
theory  for  the  isothermal  wall.  A  comparison  of  Sutherland's  law 
and  a  low  temperature  viscosity  calculation  showed  little  difference 
for  the  high  speed  flow.  Theory  based  on  reference  temperature 
supported  this  result.  The  program  was  modified  for  wall  temperature 
step  over  the  insulating  article  to  model  heat  transfer  coefficient 
dependence  on  wall  temperature  variation.  The  step  size  was  varied 
by  doubling  and  quadrupling  the  initial  step  size,  and  by  doubling 
the  initial  wall  temperature.  The  magnitude  of  the  heat  transfer 
coefficient  decreased,  and  response  slov;ed,  with  increasing  step 
size.  The  step  location  was  also  varied.  Analysis  by  v,^all  tempera¬ 
ture  superposition  yielded  results  similar  to  program  results  v/ith 
the  initial  wall  temperature  step.  The  nonisothermal  wall  effect 
appears  to  account  for  most  of  the  discrepancy  in  the  original  wind 
tunnel  test  results. 


A 


tECUfllTV  CLAStiriCATlON  Of  THIS  f  AOCfVAMi  Om 


